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(54) Title: METHOD FOR CONTROLLING CHEMICAL SYNTHESIS PROCESSES 
(54) TItre: METHODE DE REGULATION DE PROCEDES DE SYNTHESE DE PRODUTTS CHIMIQXJES 
(57) Abstract 

A method for controlling a process for synthesising at least one chemical in 
a plant including at least one reactor (R) that may be considered as a perfectly 
mbced reactor, wherein control variables (GQ enable the process to be controlled 
so tiiat a number of variables known as controlled variables (GR) and related to the 
properties of the product and/or the operation of the process match corresponding 
set values (Cgr). Said method comprises the steps of (a) inputting set values for the 
controlled variables (Cqr); (b) using a prediction unit (OP) to calculate predictions 
of the controlled variables (Pqr) on the basis of measurements of the process control 
variables (Mqc); (c) using a control unit (OC) to calculate set values for the process control variables (Coc) on the basis of the set values 
(Cgr) and predictions (Pgr) of the controlled variables; and (d) transmitting set values for the process control variables (Cgc) to process 
control actuators. The prediction unit (OP) is based on a mathematical process model known as a direct model (M), and designed in such 
a way diat the mass Mxr of at least one component (X) in the reactor (R) is predicted by means of the equation Mxr - LAG (FxRin . rx, 
rx). wherein FxRin is the mass flow rate of component X at the inlet of the reactor R, rx is the dwell time of x in the reactor, and die 
function y - LAG (u, r) is the solution of the differential equation r . dy/dt + y. calculated with the instantaneous value of u and r as 
well as with the last calculated value of y. 




(57) Abr^^ 

Une mdthode de regulation d'un procfdd de synthase d'au moins un produit chimiquc dans un appareiUage comprenant au moins un 
rdacteur (R) pouvant 6tre assimil6 ^ un rdacteur parfaitement m61ang6, dans laquelle des grandeurs de commande (GC) peimettent d'agir sur 
le dferoulement du proc6d6 en vue que des ou plusieurs grandeurs litfes aux propri6tes du produit et/ou au deroulcment du procddd. appeltes 
grandeurs i^gl^es (GR), soient 6gales & des consignes correspondantes (Cgr), ladite mdthode comprenant les Stapes suivantes: (a) entree 
de consignes concemanl les grandeurs rfigldes (Cgr); (b) calcul, au moyen d'un organe de prediction (OP), de predictions des grandeurs 
i^giecs (Pgr)» sur la base de mesurcs des grandeurs de commande du proc^d^ (Mgc); (c) utilisation d*un organe de contr61e (OC) pour 
calculer des consignes des grandeurs de commande du procede (Coc), sur base des consignes (Cgr) et des predictions (Pgr) des grandeurs 
r^giecs; (d) transmission des consignes des grandeurs de commande du procede (Cgc) ^ des actionneurs afin d'agk sur le deroulement du 
procede; dans laquelle Torgane de prediction (OP) est base sur un module matii^matique du proc^de, appeie module direct (M), et est congu 
de telle manifere qu'on predit la masse Mxr d'au moms un constituant (X) dans le reacteur (R) par requation: Mxr - LAG (FxRin . rx, 
rx), dans laquelle FxRin est le debit massique du constitaant X entrant dans le reacteur R; rx est le temps de sejour de x dans le reacteur, 
la fonction y - LAG (u, r) est la solution de ^equation differentielle r . dy/dt + y. calcuiee avec la valeur instantanee de u et de r. ainsi 
qu'avec la demi^re valeur de y calcuieee. 
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M6thode de regulation de proc6des de svnthese de produits chimiaues 



La presente invention concerne une methode de regulation de 
procedes de synthase de produits chimiques. EUe concerne egalement un 
dispositif de regulation pour la mise en oeuvre de cette methode, ainsi qu'un 
proc&l^ de synthese, en particulier de polymere, regule par ceue methode. 

S Dans un precede de synthese de produits chimiques conduit 

conventionnellement, des regulateurs de type PID (propoitionnel - integral - 
differential) sont utilises pour reguler individuellement un nombre plus ou 
moins important de grandeurs (temperatures, debits, pressions, .,.) ayant 
une influence sur le deroulement de la synthese. En d'autres termes, pour 

10 chaqiie temperature, debit ou pression a reguler, on mesure en continu (ou 
par interraittence) sa valeur effective, et un regulateur PID compare cette 
valeur effective k une valeur de consigne et agit sur la grandeur a 
commander de fagon k rdduire, le cas echdant, la difference entre la valeur 
de consigne et la valeur mesuree. 

15 Vu la complexity de la plupart des procedes industries de synthase 

chimique, les valeurs de consignes des differents regulateurs doivent 
aujourd'hui encore £tre ajustees empiriquement pour obtenir finalement les 
propriet^s souhaitees du produit synthdtise. On utilise a cette fin des 
recettes, qui foumissent des combinaisons de parametres determinees 

20 empiriquement pour obtenir, en regime stable, les propriety souhaitees du 
produit synth^tise. 

De ces recettes pourraient etre deduites, k Taide d'outils statistiques 
plus ou moins sophistiques, des relations empiriques entre les grandeurs 
reglees et les proprietes du produit synthetise. II est cependant evident que 

25 ces relations empiriques ne peuvent guere tenir compte des multiples 

interdependances existant entre les differentes grandeurs reglees separement, 
pas plus que des perturbations inconnues telles que les teneurs en impuretes 
des matieres premieres. 

n est aussi evident qu'une regulation classique a boucle fermee, uti- 

30 lisant comme corrections de retour des mesures de proprietes essentielles du 
produit synthetise, est difficilement applicable a la plupart des procedes de 
synthese. En effet, les temps morts intervenant soit dans le proc&le, soit 



dans les mesures ou analyses utUisees comme corrections de rctour, sont 
trop eleves, et les interdependances entre les diff^rcntes grandeurs regissant 
le precede sont trop complexes. 

II est connu depuis longtemps que les procedes de synthese, 
notamment les procedes de synthese en continu de polymeres (procedes de 
polymerisation), equipes de regulations avec ajustage empirique des valeurs 
de consignes, ont des inconvenients majeurs, qui peuvent etre resumes 
comme suit : 

- les demarrages du procede de synthase prennent beaucoup de temps et 
g^nfenent des quantites importantes de pnoduit "hors normes" ; 

- les transitions de qualitd sont lentes« ce qui entraine egalement la produc- 
tion de quantites importantes de produits de transition "hors normes" ; 

- Tallure de marche du procede, c'est-a-dire la masse de produit(s) 
synthdtisee par unite de temps, peut difficilement etre modifiee sans 
alterer les proprietes de ce ou ces produits ; 

- la Constance des propri6t& essentielles du ou des produits synthetises 
laisse souvent k desirer, mgme en regime stable. 

Afin d'6viter un ajustage empirique des valeurs de consigne, il a 6i6 
propose dans la litterature specialisee d'utiliser des methodes de regulation 
de procedes de synthese qui ont recours k des equations caracteristiques 
modelisant le procede de synthese pour relier certaines proprietes du ou des 
produits synthetises aux conditions de marche du ou des reacteurs lors de la 
synthase. Cependant, afin de limiter la complexite de ces equations caracte- 
ristiques, on considdrait jusqu'^ present qu'on etait obligd en pratique soit 
de consid6rcr exclusivement le cas statique (regime stable), soit de se 
limiter a une modelisation empirique tres simplifife de la dynamique du 
procede. L'utilisation d'un modfele statique est limitee au contiole d'un 
regime de production assez stable. 

Dans le cas d'une modelisation empirique, les equations 
caracteristiques ne sont valables que pour une plage de vaiidite etroite (a 
proximity du point oii la modelisation a et6 effectuee). Dans les deux cas, 
les phases de ddmarrage et les phases de transition sont mal maitrisees. II 
est certes envisageable de "couvrir" une plage plus large de conditions de 
fonctionnement en procedant a une telle modelisation locale en plusieurs 
points distincts de Tespace des parametres de fonctionnement, mais une telle 
approche devient prohibitive des que Ton tente de reguler plusieurs gran- 
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deurs en agissant sur de mult^Ies panimetres. 

II serait par consequent souhaitabie de disposer d'une methode et 
d'un dispositif de regulation simples, mieux adapt& aux specificites de la 
dynamique de procedes de synthese de produits chimiques. 

A cette fin, la presente invention conceme une methode de regulation 
d'un precede de synthase d'au moins un produit chimique dans un 
appareillage comprenant au moins un reacteur (R) pouvant etre assimile a 
un reacteur parfaitement melange, dans laquelle une ou plusieurs grandeurs 
de commande (GC) permettent d'agir sur le deroulemeni du procede en vue 
qu'une ou plusieurs grandeurs li&s aux proprietes du produit et/ou au 
d^roulement du procdde, appelees grandeurs reglees (GR), soient 6gales k 
des consignes correspondantes (Cqj^) (ou du moins aussi proches que 
possible de celles-ci), ladite methode comprenant les etapes suivantes : 

(a) entree de consignes concemant les grandeurs reglees (Cqj^); 

(b) calcul, au moyen d'un organe de prediction (OP), de predictions des 
grandeurs reglees (Pqr), sur la base de mesures des grandeurs de 
commande du pn)c6d6 (Mqq); 

(c) utilisation d'un organe de contr61e (OC) pour calculer des consignes 
des grandeurs de commande du procedd (Cq^), sur base des 
consignes (Cqj^) et des predictions (Pqr) des grandeurs reglees; 

(d) transmission des consignes des grandeurs de commande du precede 
(Cqq) k des actionneurs, ou a des organes de reglage controlant des 
actionneurs, afm d'agir sur le d^roulement du procede; 

dans laquelle Torgane de prediction (OP) est base sur un modfele 
mathematique du procede, appele modele direct (M), et est congu de telle 
manibre qu'on prddit la masse M^r d'au moins un constituant (X) dans le 
reacteur (R) par Tequation : 

^XR = LAG (FxRin • ' > 

dans laquelle : 

- PxRin ^^^^^ massique du constituant X entrant dans le reacteur R ; 

- Ty^ est le temps de sejour de x dans le r&icteur (constante de temps), qui 
vaut 

^X = '^XR ^ F^dis) 

dans lequel : 

- designe la demiere valeur calcul& de la masse du 

constituant X presente dans le reacteur R ; 
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- E Fxdis designe la somme de tous les debits riiassiques Fxdis avec 

lesquels le constituant X disparait du reacteur R, notamment par 

reaction et/ou par sortie du reacteur ; 

- la fonction y = LAG (u , r) est la solution de Tequation differcntielle 

dy 

u = T • + y 

5 calculee avec la valeur instantanee de u et de t, ainsi qu'avec la deraiere 
valeur de y calculee. 

L'avantage de cette methode est que Tequation differcntielle 
susmentionnee est resolue par un calcul algebrique simple, par exemple par 
la formule suivante (T designant I'intervalle de temps, generalement petit 
10 par rapport a r, separant les calculs successifs) ou par une formule 
equivalente a celle-ci : 

y(t-T) + u(t)'—L^ 
y(t) = lllL 



1 + 



T(t) 

Dans le cas oi^ les masses de plusieurs constituants sont evalu^es 
comme expose ci-dessus, la methode de Tinvention est particulierement 
avantageuse dans la mesure ou ces masses peuvent Stre calculees 

15 sequentiellement par de tels calculs algebriques simples, recalcul& 
frequemment (T < < r en general). Par opposition, les methodes 
traditionnelles necessitent la resolution simultanee d'un systeme d'equations 
differentielles, ce qui exige gen^raleihent une grande puissance de calcul et 
des algorithmes de calcul numdrique (integration) sophistiques ; il en resulte 

20 que la duree de chaque iteration de calcul est elevfe, et par consequent une 
regulation de ce type reagit mal k des variations rapides. 

Le precede de synthase regule peut servir a la synthase d'un 
compose monomerique ou polymerique ; de tr^s bons resultats ont 6te 
. obtenus pour la regulation de precedes de polymerisation. Le precede 

25 s*etend aussi au cas ou plusieurs produits utiles sont synthetises 

simultanement dans le meme prooede. . Le precede peut etrc continu ou 
discontinu (batch) ; la methode de regulation de Tinvention donne 
d'excellents resultats dans le cas de proc6des continus. Le precede de 
synthese regule peut eventuellement ne constituer qu'une partie d'un 

30 precede plus vaste, dent les autres parties sont reguldes d'une maniere 

differenie, ou ne sont pas reguiees. Pour que la methode de regulation de 



wo 96/41822 PCT/EP96/02408 



-5 - 

rinvention soit applicable, il faut qu'au moins un reacteur puisse etre 
assimile a un reacteur parfaitement melange, c'est-i-dire k un r&cteur dans 
lequel les differentes grandeurs (temperature, concentrations des constituants 
presents, etc.) sont quasiment identiques en chaque point. D'autres 
5 neacteurs eventuels peuvent etre du type piston ("plug-flow") : on les 
modelise mathematiquement par des temps morts. La methode est 
egalement applicable a un procede se deroulant dans plusieurs reacteurs 
disposes en serie et/ou en parallele, pouvant produire des produits aux 
propriety identiques ou diffiSrcntes. 

10 Par "constituants", on entend designer Tensemble des substances 

presentes dans le reacteur et destinees k participer k la synthese ou i la 
permettre : non seulement les reactifs de depart ainsi que le ou les produits 
synthetis6s, mais egalement la ou les eventuelles substances ne subissant 
aucune transformation, telles que solvants, catalyseurs, etc. 

15 Outre un ou plusieurs reacteurs, Tinstallation dans laquelle se 

deroule le procede peut eventuellement comprendre d'autres dispositifs 
classiques tels que ddtendeurs, strippers, condenseurs, secheurs, colonnes de 
distillation, etc. G6n6ralement, ces dispositifs auxiliaires peuvent d'ailleurs 
etre Egalement considdr6s comme des reacteurs (parfaitement melanges ou 

20 de type piston), meme si aucune reaction chimique ne s'y deroule. 

Dans le cas d'un proc6d^ de polymerisation, les "grandeurs liees aux 
proprietes du produit" peuvent par exemple etre choisies parmi la masse 
moleculaire, Tindice de fluidite en fondu, la masse volumique standard, la 
teneur en comonom&re lorsqu*un comonom&re est present, etc. 

25 Des exemples de "grandeurs lides au proced^" sont notamment la 

temperature et la pression regnant dans le reacteur, Failure de marche du 
proc6de, les concentrations des difFerents rdactifs dans le reacteur, etc. 
L'allure de marche designe la masse de produit synthetisee par unit6 de 
temps, qui n'est cependant pas forcement egale au debit de produit 

30 synth6tise sortant du reacteur : ainsi, a titre d'exemple, dans les phases de 
demarrage notamment, le debit massique de produit synthetise sortant du 
reacteur est gen^ralemem uts faible, voire nul, bien que la synthase ait 
commence, c*est-k-dire que ce debit sortant est alors inf(£rieur k Tallure de 
marche. En regime stable, par contre, on peut assimiler Failure de marche 

35 k la masse de produit synthetise par unit6 de temps. 

Des exemples de "grandeurs de commande!' sont les debits de 
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reactifs entrant dans le reacteur, la puissance foumie aux dispositifs de 
chauffage, etc. Ce sont ces grandeurs qui permettent d*agir sur le 
ddroulement du proced^ ainsi que, generalement, sur les propriet^s du 
produit synth6tis6. 

5 La ou les consignes de la ou des grandeurs de commande sont 

transmises directement ou indirectement k des actionneurs classiques tels 
que notamment des vannes, des elements chauffants, etc. "Indirectement" 
signifie que les grandeurs de commande peuvent etre transmises par 
rintermediaire d'un ou plusieurs organes de regulation (contrdlant 

10 generalement une seule variable, par ex. des regulateurs PID) contrdlant le 
ou les actionneurs (regulation "locale"). 

Sur le plan materiel, Torgane de prediction et I'organe de controle 
sont generalement des dispositifs de calcul classiques permettant d'effectuer 
des calculs en fonction de leur cablage ou de leur programmation ; il peut 

15 notamment s'agtr d'ordinateurs ou de systemes numeriques de controle- 

commande (SNCC). Un seul dispositif peut avantageusement combiner les 
fonctions de prediction et de contr6Ie. Le ou les dispositifs de calcul 
utilises sont de preference du type numerique, et foumissent periodiquement 
(par intermittence) les r6sultats de leurs calculs. Les intervalles de temps 

20 s^parant la foumiture de ces resultais peuvent varier dans le temps, et 

peuvent egalement differer selon le resultat conceme : il est clair que des 
grandeurs k variation rapide doivent etre recalculees plus frequemment que 
des grandeurs k variation lente. Des registres li d6calage peuvent etre 
utilises pour simuler materiellement les temps morts. 

25 L*organe de prediction est base sur un modele mathematique direct 

du procede (M), dans lequel le reacteur (R) est assimile k un reacteur 
parfaitement melange ; un ou plusieurs delais purs (temps morts) peuvent 
eventuellement etre pris en compte pour representer d'eventuels reacteurs de 
type piston, d*eventuels delais de transport ou delais d'obtention de resultats 

30 de mesures, .... 

L*organe de contrdle est de preference base sur Finverse du modele 
direct utilise dans Torgane de prediction (modele inverse). 

En general, la somme E Fxdis de tons les debits massiques (Fxdis) 
avec lesquels le constituant X disparait du reacteur R comprend deux 

35 termes : 

FRj^ qui designe le debit massique suivant lequel X est consomme par 
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une ou plusieurs reactions chimiques dventuelles ; 
* ^Xout designe T^ventuel debit massique de X sonant du reacteur 
par soutirage en cours de reaction, dans le cas (habituel) ou X n'est pas 
entierement consomme par reaction dans ce reacteur ; ou encore, par 
5 exemple, par evaporation dans le cas d'un reacteur ouveit. 

L'interet de la methode est que les termes Fxdis sont en general 
proportionnels k ' exemple, on a en general 
^Xout ^^XR / 

(tjj designant le temps de sejour du reacteur R) 
10 et Frx = Rx*Mxr 

(Rx designant la reactivite de X dans le reacteur R). 

Dans ce cas, Texpression donnant se simplifie et devient : 

rx = 1 / (Rx + I/tr) 

Cette expression est independante de Mjqj, ce qui constitue une 
15 simplification extremement interessante. 

Un autre avantage de la mdhode reside dans le calcul periodique du 

temps de sejour r^. En effet, r^resente bien la dynamique du 

constituant consider^ dans le reacteur. Ceci peimet notamment de suivre 

revolution de ce parametre, ce qui est important pour la comprehension de 
20 la dynamique du procede, et par consequent pour sa regulation. Au 

contraire, les methodes empiriques de type "boite noire" ne pennettent pas 

d'acceder a ce parametre. 

Avantageusement, le calcul des predictions des grandeurs reglees 

(Pqr) P®^t additionnellemeni tenir compte d*une ou plusieurs mesures de 
25 grandeurs reglees (Mqr), de grandeurs de commande (Mqq) et/ou d'autres 

grandeurs liees au deroulement du proc^dd (M^). 

De meme, avantageusement, le calcul des consignes des grandeurs 

de commande du procede (Cq^) peui additionnellement tenir compte d'une 

ou plusieurs mesures de grandeurs reglees (Mqr), de grandeurs de 
30 commande (Mqq) et/ou d'autres grandeurs liees au deroulement du procede 

(M^), identiques ou diffSrentes de celles eventueUement prises en compte 

pour le calcul des predictions des grandeurs reglees (Pqr)' 

Toutes les mesures dont il est question dans la presente description 

ne sont pas forcement des mesures directes, dans le sens ou Tune ou 
35 plusieurs d'entre elles peuvent eventueUement etre des mesures 

inferentielles, c'est-a-dire des valeurs obtenues par calcul a partir d'une ou 
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piusieurs autres mesures directes. Ainsi, par exemple, Talluie de maiche 
de certains precedes de synthese exothermiques ne peut pas se mesurer 
directement, mais on peut en obtenir une mesure inferentielle par calcul, par 
exemple a partir de mesures (directes) du debit et des temperatures d'entree 
et de sortie du fluide de refroidissement. 

Dans le cas particulier de procedes de polymerisation, la ou les 
proprietes du polymere intervenant dans la regulation sont de prefi^rence 
choisies parmi la masse volumique specifique (MVS) du polymere, les 
proprietes rh^ologiques du polymere k T^tat fondu, et sa teneur en 
comonom&re. En particulier, la ou les proprietes rheologiques intervenant 
dans la methode de regulation sont avantageusement Tindice de fluidite en 
fondu (melt index) du polymere et/ou une mesure de viscosite. 

Avantageusement, on evalue une ou plusieurs propriety du polym&re 
en utilisant une technique choisie parmi la spectroscopie infrarouge proche 
(NIR), la spectroscopie infrarouge k transformee de Fourier (FUR) et la 
resonance magn6tique nucleaire (NMR). 

En particulier, on peut avantageusement ^valuer une ou plusieurs 
proprietes du polymere en s^liquant une relation de correlation preetablie 
aux resultats de mesures realisees par spectroscopie infrarouge proche (NIR) 
a plusieurs longueurs d'onde predeterminees en fonction de la nature du 
polymere, choisies entre 0,8 et 2,6 mm. 

Davantage de details concemant la realisation de telles mesures dans 
le cadre de ia regulation de procedes de polymerisation peuvent etre trouves 
dans la demande de brevet EP 328826 (US S1SS184). 

En vue de tenir compte d'ecaits eventuels entre les mesures et les 
predictions des grandeurs regiees, il peut etre utile de recourir k une 
correction. 

Un premier type de correction consiste en ce que la consigne d'au 
moins une grandeur reglee (Cqj^) soit corrigee sur la base de l^ecart 
(avantageusement filtre) entre la mesure (Mqj^) et la prediction (Pqr) de 
cette grandeur reglee, de fagon k ce que la regulation soit efficace (Mq^^ = 
Cqj^) meme en presence d'une erreur dans la prediction de cette grandeur 
reglee. Cette technique est couramment designee par Texpression "internal 
model control" (IMC). 

Un second type de correction consiste en ce que le modele (M) du 
procede soit periodiquement adapte sur la base de Tecart (avantageusement 
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filtre) entre les predictions (Pqr) mesures (Mqj^) des grandeurs 
rcgl6es, de fagon a ce que, ici egalement, le modele du precede foumisse 
des predictions des grandeurs r6gl&s (Pgr) aussi proches que possible 
(idealement, egales) des mesures de ces grandeurs (Mq^^), ce qui est 
5 indispensable a une regulation efficace. 

L'adaptation consiste k recalibrer le modele, c'est-a-dire a en 
recalculer un ou plusieurs parametres ; normalement, le nombre de 
parametrcs recalcules ne depasse pas le nombre des grandeurs reglees pour 
lesquelles on dispose k la fois d'une prediction et d'une mesure. Une 

10 resynchronisation (recalage dans le temps) de ces mesures est sou vent 
souhaitable, surtout lorsqu'il s'agit de mesures de proprietes du produit 
synthetise dont la duree d'obtention est longue. Ce second type de 
correction est plus avantageux dans la mesure b& il permet d^adapter le 
modele egalement sur le plan de sa dynamique. 

15 L'adaptation conceme non seulement le modele direct du procede 

(organe de prediction), mais egalement le modMe inverse (organe de 
contr&le), 

Selon une variante avantageuse, les mesures (Mqj^) des grandeurs 
r6gl6e$ interviennent uniquement dans l'adaptation ^ventueUe du modele du 
20 procede, et n' interviennent pas directement dans le calcul des consignes des 
grandeurs de commande du procede (Cqq), 

C'est-a-dire que les mesures des grandeurs reglees n' interviennent 
pas dans la regulation proprement dite : Tavantage en est que la qualite de 
la regulation n'est ainsi pas affectde par la lenteur 6ventuelle de revaluation 
25 des proprietes du produit. 

Un autre aspect de I'invention conceme une methode de regulation 
telle que decrite ci-dessus, appliquee a un procede de polymerisation, 
comprenant Tune ou plusieurs des Stapes supplementaires suivantes : 

calcul d'une consigne de temperature dans le reacteur en fonction d'une 
30 ou plusieurs consignes de proprietes du produit ; et transmission de 

cette consigne de temperature a un ou plusieurs actionneurs permettant 
de modifier la temperature dans le ri6acteur (6ventuellement - 
indirectement, c'est-a-dirc par I'interm&iiaire d'un ou plusieurs organes 
de regulation, par ex. des regulateurs PID, contr&lant le ou les 
35 actionneurs) ; 

calcul d'un bilan thermique pour le reacteur, sur base notamment de 
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mesures de temperature; utilisation de ce bilan thermique de fa^on k 
determiner la quantite de polymere synthetise par unite de temps (allure 
de marche) et/ou la productivity du catalyseur et/ou la concentration 
d'au moins un reactif dans le reacteur ; 
5 - calcul de la quantite de chaleur produite par la polymerisation, grace a 
un calcul de la quantite du ou des reactifs qui polymerisent ; 
-determination par ce biais de la quantite la chaleur qu'il faut ajouter ou 
evacuer pour maintenir la temperature du reacteur ; utilisation du 
resultat dudit calcul (par exemple par feed-forward) pour ameliorer la 
10 regulation de temperature, de fagon k respecter au mieux la consigne de 

temperature, notamment en cas de modifications de Failure de marche. 
Ces variantes sont basees sur le lien qui existe entre la quantite du 
ou des reactifs participant a la reaction et la quantite de chaleur produite ou 
absorbee par la reaction. 
15 Selon une variante avantageuse, la propriete Pxj^ d'un constituant 

"x" dans le reacteur R, assimile k un reacteur parfaitement melange, est 
calculee comme suit : 

Pxr = LAG (PxjN , Mxr / Fxijj ) 
ou "Px" est une propriete d'un constituant "x", repondant de manibre 
20 substantielle a la loi de melange lineaire Pxi^2 ^1 • ^1 ^2 • 

Px2 , 

wj et W2 etant les proportions massiques de deux fractions 1 et 2 de 

propriete Px^ et Px2 que Ton melange (wj + W2 = 1) ; 
Pxj ^2 ^ propriety de x a sa sortie du reacteur, aprfes melange ; 
25 ^IN ^ propriete du constituant "x" h son entree dans le reacteur R ; 
Mx]^ est la masse du constituant x dans le reacteur R ; 
Fxjjj est le debit massique du constituant x entrant dans le reacteur R. 

Une transformation mathematique permet parfois de rendre lineaires 
(additives) certaines grandeurs qui ne le sont pas : par exemple, Tindice de 
30 fluidite en fondu d'un polymere (melt index) ne rdpond pas a une loi de 
melange lineaire, mais bien son logarithme ; le calcul de Px| ^2 
susmentionne s'effecue des lors sur le logarithme de ce parambtre. 

Selon une autre variante avantageuse, la methode de regulation de 
rinvention comprend les etapes suivantes : 
35 - entree de consignes relatives k une ou plusieurs proprietes du produit a 
synthdtiser, dans un algorithme maitre ; 
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- entree de la consigne d 'allure de marche du precede, dans un algorithme 
esclave; 

- calcul des consignes de concentration des constituants dans le rfacteur k 
I'aide de ralgorithme maitre, en fonction notamment des consignes et des 

5 . mesures de proprietes du produit ainsi que de mesures ou de pr&iictions 
des concentrations des differents constituants dans le r&cteur ; 

- transmission des consignes de concentration calculees par ralgorithme 
maitre comme grandeurs d'entree a ralgorithme esclave; 

- calcul de consignes de debit des constituants entrant dans le reacteur, k 
10 I'aide de I'algorithme esclave, en fonction notamment de la consigne 

d'allure de marche du proc6de, des consignes de concentration, et de 
mesures de debit des constituants entrant dans le reacteur, et 

- transmission des consignes de debit calculees a Taide de Talgorithme 
esclave a un ou plusieurs actionneurs (^ventuellement indiiectement, 

15 c'est-i-dire par Tinterm^diaire d'un ou plusieurs organes de regulation, 

par ex. des rcgulateurs PID, controlant le ou les actionneurs) afin de 
r^gler les debits des constituants entrant dans le reacteur, 
dans laquelle ralgorithme msdtre et/ou I'algorithme esclave sont mis en 
oeuvre comme decrit ci-dessus, c'est-^-dire en utilisant la fonction LAG 
20 pour calculer la masse d'au moins un constituant dans le reacteur. 

Les algorithmes maitre et esclave sont egalement mis en oeuvre au 
moyen d'un ou plusieurs dispositifs de calcul classiques. Selon une variante 
avantageuse, Tensemble des calculs (prediction, controle, etc.) de ces deux 
algorithmes sont effectues par un meme dispositif de calcul. 
25 Avantageusement, des mesures de temperature (par ex. temperature 

dans le reacteur, et/ou temperature d'entree et/ou de sortie d'un eventuel 
fluide de refroidissement) interviennent comme grandeurs d'entree 
supplementaires dans Torgane de prediction et/ou de controle. 

De preference, ralgorithme esclave prend Egalement en compte des 
30 mesures de la composition de constituants presents dans ou sortant du 
reacteur. 

De maniere avantageuse, la m^thode de regulation comprend en 
outre une 6tape de calcul a I'aide de I'algorithme esclave, en fonction des 
mesures de debits, de predictions de concentrations transmises k 
35 ralgorithme maitre pour calculer des predictions de proprietes servant 

comme grandeurs d'entree supplementaires dans le calcul des consignes de 
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concentrations. 

Les algorithmes maitre et esclave forment une regulation du type 
cascade. H est particulibrement avantageux que raigorithme maitre et/ou 
ralgorithme esclave soient adaptatifs, c'est-a-dire que certains de leurs 

5 parametres soient periodiquement recalcules intervalles reguliers ou 
irreguUers). Une telle adaptation permet notamment de garantir que le 
modele mathematique decrit le plus fidelement possible le precede dans son 
etat courant, meme en cas de modification de certaines conditions 
operatoines (temperature, pression, allure de marche, etc.) et en cas de 

10 perturbations (empoisonnement du catalyseur, ...). 

L'algorithme maitre r^se le contrdle des propriet6s du produit & 
I'aide d'un modele base sur des equations caracteristiques reliant les 
proprietes du produit aux concentrations des differents constituants dans le 
reacteur, ainsi qu'eventuellement a la temperature regnant dans le reacteur. 

15 L'algorithme esclave regule les concentrations d'un ou plusieurs constituants 
en agissant sur les debits d'alimentation d'un ou plusieurs constituants 
6ventueUement differents. 

L'intergt de cette cascade "maitre-esclave" reside dans le fait que le 
module maitre determine avec precision les concentrations de constituants 

20 necessaires a Tobtention des proprietes voulues pour le produit synthetise, et 
que le modele esclave assure le respect des valeurs imposees par le maitre. 
Pilote par le maitre, Tesclave est par consequent capable : 
- d'amener rapidement les concentrations aux valeurs desirees par le maitre 
et de les y maintenir, 

25 -de controler efficacement Failure de marche du precede sans peiturber 
les concentrations. 

Cette cascade maitre-esclave est paiticulierement efficace car aussi 
bien le maitre que T esclave prennent en compte la dynamique du precede 
grace k Tutilisation de la fonction LAG dans les calculs. 

30 L'algorithme esclave peut en outre etre con^u de maniere a foumir a 

ralgorithme maitre des predictions de concentrations fiables. De ces 
predictions ou de mesures de concentrations, Talgorithme maitre deduit des 
predictions de proprietes fiables du produit en train d'fetre synthetis6 dans le 
reacteur. En comparant ces predictions de propriety aux consignes de 

35 proprietes, ralgorithme maitre peut. le cas echeant, intervenir et corriger 
les consignes de concentrations. Cette correction est possible avant meme 
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qu'un ecart entre une variable et sa consigne ne se presente. La prise en 
consideration de predictions de propriet^s obtenues a partir de predictions 
ou de mesures de concentrations permet de reduire considerablement les 
fluctuations temporelles des proprietes du produit synthetise et il en resulte 
une meilleure Constance de la qualite du produit. 

Si les proprietes du produit a synthetiser sont dependantes de la 
temperature dans le ou les rfiacteurs, il est preferable de prevoir une 
regulation de la temperature par Talgorithme esclave. Ce dernier dresse le 
bilan thermique de chaque reacteur et determine, au moyen du calcul de 
I'aUure de marche, la quantite de chaleur qu'il faut ajouter ou 6vacuer pour 
respecter des consignes de temperature calcul^es par Talgorithme mdtre. 
De ces resultats il deduit des consignes d'entree pour des organes de 
regulation thermique de TappareiUage de synthese. II sera apprecie que 
cette fa?on de proceder permet d'intervenir sur les organes de regulation 
thermique de Tappareillage de synthese avant meme que la temperature ne 
change. Des mesures de temperatures interviennent par ailleurs 
avantageusement comme grandeurs d'entree suppiementaires dans 
Talgorithme esclave. 

L'algorithme maitre comprend avantageusement la structure 
suivante : 

- un organe de prediction, base sur un modele direct du procede permettant 
de foumir une prediction des proprietes du produit synthetise en fonction de 
mesures et/ou de predictions des concentrations des constituants ; 

- un organe d'adaptation comparant les predictions de proprietes calcuiees 
par Torgane de prediction a des valeurs effectivement mesurees sur le 
produit synthetise et deduisant de cette comparaison des parametres 
d*adaptation, lesdits parametres d'adaptation intervenant comme grandeurs 
d'entree suppiementaires dans ledit organe de prediction de Talgorithme 
maitre ; et 

- un organe de controle, base sur un modele inverse du precede, pour 
calculer, en fonction des consignes et des predictions de proprietes du 
produit k synthetiser, des consignes de concentrations pour Talgorithme 
esclave, lesdits parametres d'adaptation intervenant egalement comme 
grandeurs d'entree suppiementaires dans ledit organe de controle. 

L'algorithme esclave comprend avantageusement la structure 
suivante : 
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- un organe de prediction, base sur un modele direct du procede permettant 
de foumir une prediction des concentrations d*un ou plusieurs des 
cx)nstituants sur base d'un bilan de matiere dans le reacteur ; 

- un organe d'adaptation, comparant les pr6dictions de concentrations 

S calculees par le modele direct a des mesures de concentration et deduisant 

de cette comparaison des parametres d'adaptation, lesdits parametres 
d'adaptation intervenant comme grandeurs d'entree supplementaires dans 
ledit organe de prediction de Talgorithme esclave ; et 

- un organe de controle, bas6 sur un modele inverse du procede, pour 

10 calculer, en fonction de la consigne d'allure de marche, des consignes de 
concentration calcul6es par roi;gane de contrdle de Talgorithme maitre et 
des predictions de concentration calculees par Torgane de prediction de 
Talgorithme esclave, les consignes pour les debits entrant dans le 
reacteur, lesdits parametres d'ad^^tation intervenant comme grandeurs 

IS d'entree supplementaires dans ledit oigane de controle de Talgoritlune 

esclave. 

La dynamique du proc6dd est avantageusement d^rite et calcul^e au 

moyen de fonctions du type y = LAG (u , r) 

cette fonction etant la solution de Tequation differentielle 

u = t • + y 
at 

20 oil les arguments u et r varient dans le temps. L* utilisation de cette 

fonction conformement aux theoremes 1 et 2 exposes ci-dessous permet de 
resoudre de maniere sequentielle les bilans de matiere utilises par 
Talgorithme esclave et de decrire la cinetique du procede par des equations 
caracteristiques simples dans I'algorithme maitre. La fonction LAG permet 

25 encore de reduire considerablement le volume des calculs n^cessaires et 
rend par consequent inutile Tutilisation d'ordinateurs rapides et puissants. 
De plus, cette fonction permet d'etablir de fagon particulierement simple les 
modeles directs et inverses du procede ou de certaines de ses parties. 
Les principales qualites de la regulation proposee peuvent etre 

30 resumees comme suit : 

- anticipation : la regulation commence a corriger les perturbations 
mesurees avant meme que leur effet ne se soit manifeste sur les mesures 
des proprietes (utilisation de predictions de concentrations, de predictions 
de proprietes et de predictions de temperatures dans les algorithmes); 
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- precision meme en presence de perturbations : le modele direct et le 
modfele inverse sont recalibres en permanence en utilisant les mesures des 
proprietes (adaptation); 

• validity etendue : Talgorithme garde sa validity durant les transitions 
. d'allure de maiche et de qualite, ainsi que lors des demanages et des 
5 arrets (mise en equation de la dynamique du procede, utilisation de 

predictions pour les grandeurs dont les mesures impliquent des temps 
morts impoitants) ; 

- simplicite : le developpement et la mise en oeuvre sont facilit^s grace a 
une methode originate de mise en ^uation de la dynamique du procede 

10 (fonaion LAG). 

Le pTOc666 de synth^e a reguler est ainsi modelise sous une forme 
qui est generalement qualifi6e de "modele de coimaissance" ("first principle 
model"), c'est-a-dire que son modfele est 61abore a partir d'equations 
refletant le deroulement physico-chimique detaill6 du procede. Une telle 

15 approche permet d'obtenir, au moyen d'un jeu d'equations relativement 

simple sur le plan mathematique, des resultats superieurs a ceux qui auraient 
ete obtenus au moyen d'un modele empirique de type boite noire, en 
foumissant notamment des param&tres relies a des grandeurs reeUes et une 
meilleure validity k Texterieur de Tespace dMdentification (extrapolation). 

20 La plupart des modeles empiriques utilisent des equations complexes, 

sou vent d'ordre eleve si Ton desire obtenir une simulation correcte de la 
dynamique du procede, et dont les parametres (notamment ies constantes de 
temps) doivent etre identifies pour un point de fonctionnement precis; le 
module n'est valable que dans ie voisinage imm6diat de ce point de 

25 fonctionnement. Une telle approche peut difRcilement 6tre generalisee a un 
grand nombre de points de fonctionnement dans le cas d'un procede de 
synthese chimique reel faisant intervenir de nombreuses grandeurs. 

Au contraire, selon la methode de regulation de Pinvention, on 
utilise un ensemble d'equations simples, purement statiques ; la dynamique 

30 du procede est simulee par des fonctions simples (cf . la fonction LAG ci- 
dessus). Avantageusement, les temps de sejour (constantes de temps des 
Equations) peuvent etre recalcules aussi souvent qu'on le souhaite, ce qui ne 
pose aucun problfeme etant donne la simplicite des equations. On obtient en 
fm de compte un ensemble d'6quations extrSmement simples et faciles k 

35 resoudre en temps reel, meme k frequence elevee. 
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La methode de regulation pnoposee est avantageusement applicable a 
des precedes de synthese, notamment en continu, de polymeres 
(polymerisation), et notamment a la polymerisation en continu d'oleflnes 
telles que par exemple Tdthylene ou le propylene, aussi bien en phase 
S liquide qu'en phase gazeuse. 

La presente invention conceme egalement un precede de synthese 
d'un ou plusieurs produits chimiques, regule au moyen d'une methode de 
regulation conforme a Tinvention. En particulier, de ties bons resultats ont 
ete obtenus pour la regulation d'un procede de synthase en continu de 

10 polyethylene par polymdrisation d'ethylene dans au moins un reacteur, les 
reactifs comprenant Tethylene, Thydrogfene et/ou un comonomere optionnel, 
la reaction de polymerisation ayant lieu dans un solvant en presence d'un 
catalyseur et une partie du contenu du reacteur etant preleve en permanence 
ou par intermittence. Ce procede peut aussi bien se derouler en phase 

15 liquide qu'en phase gazeuse ; de preference, il se deroule en phase liquide 
(dans un solvant). 

La methode s'applique de mani^re analogue a la synthese du 
polypropylene (le monomere de depart princqial etant dans ce cas le 
propylene au lieu de rethylene), du propane pouvant egalement etre present 

20 si le procede se deroule en phase gazeuse. Pour le polypropylene, Tindice 
de fluidite en fondu est souvent designe par MET au lieu de MI. 

L' invention conceme egalement un dispositif de regulation destine a 
mettre en oeuvre la methode de regulation de Tinvention, ainsi qu'une 
installation de synthase d'un ou plusieurs produits chimiques comprenant un 

25 tel dispositif de regulation. 

De maniere plus precise, T invention conceme egalement un 
dispositif de regulation d'un procede de synthase d'un produit chimique dans 
un appareillage de synthese comprenant au moins un reacteur, ledit 
dispositif comprenant : 

30 - au moins une unite de calcul; 

- des moyens pour entrer des consignes de proprietes du produit a synthe- 
tiser dans Tunite de calcul; 

- des moyens pour entrer des consignes d'allure de marche du produit a 
synthetiser dans I'unite de calcul; 

35 - des organes de mesure de debit des flux entrant dans le reacteur; 

- des organes de mesure de la composition des flux sortant du reacteur; 
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- des organes de commande de debits (actionneurs) pour regler les debits 
des flux entrant dans le r^cteur; 

- des moyens de communication entre ladite unite de caicul, iesdits organes 
de mesures de debit et Iesdits organes de regulation; 

5 dans lequel : 

- la masse d'au moins un constituant est calculee par la fonction LAG 
comme expose ci-dessus; 

- Tunite de caicul est apte a calculer a I'aide d'un algorithme maitre, en 
fonction des consignes de propriet^s, des consignes de concentration de 

10 reactifs dans le reacteur; 

- Tunite de caicul est apte a calculer a Taide d'un algorithme esclave, en 
fonction des consignes de production et des consignes de concentration, 
de consignes de debit pour les flux entrant dans le reacteur, ces 
consignes de debit etant transmises comme consignes d'entree aux 

IS oi^anes de regulation de debits; 

- les mesures r6alis^s par les organes de mesure de debits interviennent 
conune grandeurs d'entree suppl6mentaires dans ledit algorithme esclave, 
pour permettre a ce dernier de calculer, en fonction de ces mesures de 
ddbits, des pr^ictions de concentrations, ces predictions de 

20 concentrations intervenant dans I'algorithme maitre pour calculer des 

predictions de proprietes utilisees comme grandeurs d*entree 
supplementaires dans le caicul des consignes de concentrations. 

L'invention conceme egalement un dispositif de regulation tel que 
decrit ci-dessus, dans lequel : 
25 - Tappareillage de synthase comprend en outre : 

des organes de regulation thermique apte a controler la temperature 

dans le reacteur, et 

des capteurs de temperatures; 

- r algorithme maitre est apte a calculer, en fonction des consignes de 
30 proprietes, des consignes de temperatures pour le reacteur; 

- Talgorithme esclave est apte a : 

calculer un bilan thermique pour le reacteur, 
- r6soudre ce ou ces bilans thermiques de fagon k determiner la 

chaleur qu'il faut ajouter a, respectivement retrancher de la synthase 
35 pour respecter les consignes de temperature, et 

deduire de ce ou ces bilans thermiques des consignes d'entree pour 



f 



wo 96/41822 



PCT/EP96/02408 



- 18 ' 



les organes de regulation themiique du reacteur, et 

recevoir comme grandeurs d*entr6e suppl^mentaires les mesures 

effectuees par les capteurs de temperature. 

L'invention conceme ^galement un dispositif tel que d&^t ci-dessus, 
S dans lequel des mesures effectuees par les capteurs de temperature 

interviennent comme grandeurs d'entree supplementaires dans Falgorithme 
maitre. 

L'invention conceme egalement un dispositif tel que decrit ci-dessus, 
comprenant : 

10 - au moins un analyseur npie k foumir des mesures des proprietes inter- 
venant dans Talgorithme maitre; et 

- des moyens pour entrer ces mesures de proprietds dans Tunit^ de calcul; 
ladite unite de calcul comprenant : 

- un premier organe de prediction, bas6 sur un premier modele direct du 
15 procede permettant la prediction des proprietes du produit synthetise, en 

fonction des pr^ictions de concentrations calculees par Talgorithme 
esclave; 

- un premier organe d'adaptation, comparant les predictions de propriety 
calculdes par ie premier organe de prediction aux valeurs effectivement 

20 mesurees sur le produit synthetis6 et d6duisant de cette comparaison des 

parambtres d' adaptation, intervenant comme grandeurs d'entrfe 
supplementaires dans ledit premier organe de prediction ; et 

- un premier organe de controle, base sur un premier modele inverse, pour 
calcuier, en fonction des consignes et des predictions de proprietes, des 

25 consignes de concentrations pour Talgorithme esclave, lesdits parametres 

d'adaptation intervenant egalement comme grandeurs d'entree 
supplementaires dans ledit premier organe de contr61e. 

L'invention conceme egalement un dispositif tel que d6cni ci-dessus, 
comprenant en outre : 
30 - au moins un analyseur apte a foumir des mesures de concentration des 
reactifs; et 

- des moyens pour entrer ces mesures de concentration dans Tunite de 
calcul; 

ladite unite de calcul comprenant : 
35 - un deuxibme organe de prediction* base sur un deuxieme module direct 
permettant la prediction des concentrations, en fonction du bilan de 
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mati^re dans le reacteur; 

- un deiixieme organe d'adaptation, comparant les predictions de 
concentrations calcul&s par le deuxifeme organe de prediction aux 
mesures de concentration et d&iuisant de cette comparaison d'autres 

. parametres d'ad^tation, intervenant comme grandeurs d*entr6e 
supplementaires dans ledit deuxieme organe de prediction; et 

- un deuxieme organe de controle, base sur un deuxieme module inverse, 
pour calculer, en fonction des consignes de production, des consignes de 
concentration calculees par I'algorithme maitre et des predictions de 
concentration du deuxieme organe de prediction, des consignes pour les 
debits entrant dans le reacteur, lesdits autres parametres d*adaptation 
intervenant comme grandeurs d'entiee supplementaires dans ledit 
deuxieme organe de controle. 

Un mode de realisation concret de Tinvention est illustre sur base 
d'un procede de synthese en continu de polyethylene (PE), en se referant 
aux Figures 1 a 10. Celles-ci montrent : 

Figure 1 un schema d'un circuit de fabrication de polyethylene; 

Figure 2 un schema simplifid de la structure d*une regulation 

avancee selon Tinvention; 
Figure 3 un schema de principe de la regulation avancde appliquee 

au circuit de fabrication de la Figure 1; 
Figure 4 : un schema de principe d'un algorithme de controle 

adaptatif, tel qu'utilise dans le systeme de regulation 

avancee selon la Figure 2; 
Figure 5 : un sch6ma de la structure de Talgorithme maitre dans le 

syst&me de regulation avancee selon la Figure 2; 
Figure 6 un schema de la structure de T algorithme esclave dans le 

systeme de regulation avancee selon la Figure 2; 
Figure 7 le schema general d'une m6thode de regulation conforme 

a r invention; 

Fig. 8-10 : le schema de variantes particulieres de la methode de 
I'invention. 

Sur la figure 1, on distingue tout d'abord le procede de synthese 
proprement dit (Pr), qui peut etre contr6ie en imposant au moins une 
consigne d'une grandeur de commande (CGC) (par exemple un ou plusieurs 
debits de constituants entrant dans le reacteur, une temperature, etc.) ^ un 
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actionneur adequat (vanhe, dispositif de chauffage ou de refroidissement, 
etc.). La regulation s'effectue au moyen d'un organe de controle (OC), 
base sur le module math^matique inverse du procede, et dont le rdle 
principal est de comparer la ou les consignes des grandeurs regimes (CGR) 
5 (par exemple une ou plusieurs proprietes du produit a synthetiser et/ou une 
ou plusieurs grandeurs liees au d^rouleraent du procede) a la ou aux 
predictions de ces grandeurs (PGR). La ou les predictions des grandeurs 
reglees (PGR) sont calculees par un organe de prediction (0P)» base sur le 
module mathdmatique direct du procdde, sur la base de mesures des 

10 grandeurs de commande (MGC). On notera qu'aucune mesure de 
propriete(s) du produit synthetise n*intervient dans la regulation. 

La figure 8 repr^sente une variante de la m&hode de la figure 7, 
dans laquelle le module mathematique du procede est periodiquement adapte 
par un organe d'adaptation (OA), sur la base de Tecart (avantageusement 

15 filtre, ou traits numeriquement) entre les predictions (PGR) et les mesures 
(MGR) des grandeurs rcglees. Une resynchronisation (recalage dans le 
temps) de ces mesures et de ces predictions est souvent ndcessaire, par 
exemple lorsquMl s'agit de mesures de proprietes du produit 5ynth6tis6 dont 
la durde d'obtention est longue. L' organe d'adaptation (OA) transmet le 

20 r6sultat de ses calculs, c'est-k-dire ses instructions d'adaptation, a Torgane 
de prediction (pour adapter le modele direct du procede) et k T organe de 
controle (pour adapter le modele inverse du procede). On notera que la ou 
les mesures de propriet6(s) du produit synthetise ne sont pris(es) en compte 
que dans le processus d'adaptation, qui se deroule generalement a une 

25 frequence beaucoup plus faible que le processus normal de controle. La 
lenteur ^ventuelle de ces mesures n'a done aucun effet direct sur la qualite 
de la regulation. 

La figure 9 r^resente une autre variante de 1' invention, dans 
laquelle une ou plusieurs mesures de grandeurs reglees (Mqj^) sont prises 

30 en compte par I'organe de controle (OC), et une ou plusieurs mesures de 
grandeurs reglees (Mqj^') (eventuellement differentes) sont prises en 
compte par Torgane de prediction (OP), De meme, une ou plusieurs 
mesures de grandeurs de commande (Mq^') peuvent en outre etre prises en 
compte par I'organe de contr61e (OC), 

35 II va de soi qu'on pourrait creer une autre variante de la methode de 

rinvention en combinant les variantes des figures 8 et 9, c'est-a-dire a la 
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fois en utilisant un organe d'adaptation et en prenant en compte une ou 
plusieurs mesures de grandeurs de commande dans T organe de controle, 
et/ou une ou plusieurs grandeurs de commande dans Torgane de prediction 
et/ou dans T organe de controle. 
5 Dans la figure 10, le modele mathematique du procede n'est pas 

adapte au sens propre du terme, mais on utilise Tecan (avantageusement 
filtre) entre les mesures et les predictions des grandeurs reglees pour 
corriger les consignes des grandeurs reglees (Cqj^). En Toccurrence, cette 
correction a 6t6 ici representee comme une simple difference : on soustrait a 

10 chacune des consignes de grandeurs reglees un tenne correctif calcule par 
I'organe d'adaptation OA (qui n'est en fait ici qu'un organe de correction), 
ce qui foumit des consignes corrigees Cq^\ transmises a Torgane de 
controle OC. II va de soi que dans certains cas, la correction peut 
impliquer des operations plus complexes qu'une soustraction, par exemple 

15 une division (dans ce cas-ci, on peut cependant se ramener a une 

soustraction en considerant les logarithmes des grandeurs considerees). 
Cette methode est commun^ment appel^ Internal Model Control (IMC). 

Si Ton se r^f&re h la Figure 1, qui represente sch6matiquement un 
circuit de synthase en continu de polyethylene (PE), la polymerisation de 

20 Tethylene se fait dans un r&cteur boucle 10, en suspension dans un solvant 
ad&juat, tel que par exemple Thexane. Le processus est continu, c*est-a- 
dire que les reactifs sont injectes en continu, et qu'une partie du contenu du 
reacteur 10 est continuellement soutiree. Une pompe de circulation (non 
representee) assure Thomogeneitd du contenu du reacteur 10. 

25 Les reactifs introduits dans le reacteur sont Tethylene "Et", 

I'hydrogene "Hy" et le butene "Bt" (cf, reference 11). Un catalyseur est 
6galement injecte en continu. II importe de bien maitriser les concentrations 
des reactifs dans le reacteur, car les proprietes de la resine PE sont 
principalement determinees par les rapports de concentrations Hy/Et et 

30 Bt/Et. 

La temperature de polymerisation dans le reaaeur est un parametre 
additionnel qui influence les proprietes de la resine PE. Vu que la reaction 
de polymerisation est fortement exothermique, la temperature du reacteur 
doit 6tre r^guiee en ay ant recours k un circuit de refroidissement 12. 
35 Le reacteur 10 en fonctionnemem contient done du solvant, du poly- 

mere et des reactifs n'ayant pas encore reagi et du catalyseur. Son contenu 
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est soutire en permanence k travers la conduite de soutirage 14. Ce contenu 
soutinS eiitre dans un stripper STP 16, qui s6pare le polymere PE et les 
fluides (solvant et reactifs). Ces fluides sont vaporises par injection de 
vapeur d'eau et ^vacues dans un condenseur CD 18. Dans ce dernier le 
5 solvant est de nouveau condense avant d'etre recycle. Les reactifs qui sont 
plus legers sont separes du solvant et sont aussi recycles. Un chromato- 
graphe en phase gazeuse GC (20) place a la sortie du condenseur 18 permet 
de determiner les concentrations Hy/Et, Bt/Et des reactifs. 

Le polymere dvacu6 du stripper 16 est concentre dans une 

10 centrifligeuse CFG 22, puis s6ch6 dans un secheur k lit fluide SHLF 24, 

avant d'6tre envoye vers le fmishing pour y &tre granule. Des echantiUons 
sont preleves a la sortie du secheur 24 pour y mesurer les proprietes de la 
resine : cristallinitd (mesuree par la masse voliimique specifique "MVS") et 
proprietes rheologiques (indice de fluidite en fondu ("melt index (MI)" ou 

15 "melt flow index (MFI)") et viscosite en fondu "^12", mesuree dans un 
^oulement capillaire sous un cisaillement de lOO s~^). 

La dynamique de ce processus de synthase de PE est lente et 
complexe : 

- Le reacteur boucle 10 se comporte comme un reacteur parfaitement 
20 melange. Des lors, toute modification du debit d*alimentation d'un des 

reactifs ne se traduira que progrcssivement sur la concentration de ce 
reactif dans le reacteur. En effet, le nouveau debit doit se melanger au 
volume tout entier du reacteur 10 pour I'amener a la nouvelle concen- 
tration d'equilibre. 

25 - La mesure des concentrations des r6actifs se fait par un chromatographe 
en phase gazeuse 20; il s'agit d'un appareil discontinu, qui proc^e par 
etapes successives : prelevement d'un echantillon de gaz, analyse, puis 
elaboration des resultats. II y a done un temps mort (de 5 a 15 min.) 
entre les changements de concentration et leur mesure. 

30 - Les proprietes du polymere fabrique a chaque instant dependent principa- 
lement des concentrations des reactifs. Toute modification de ces 
concentrations affecte done instantanement les proprietes du polym&re 
fabrique. Par contre, les proprietes moyennes dans le reacteur ne se 
modifient que progressivement, car le polymere fraichement produit doit 

35 se melanger au polymere deja present dans le reacteur 10 (temps de 

sejour ; ± 2 h). 
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- Lorsque le poly mere est soutire du reacteur 10, il subit a nouveau une 
serie de melanges dans les differents appareils (STP, CFG et SHU) 
•destines a le s6cher (temps de sejoiir : ± 2 h). Des ecbantillons de 
polymere sont alors pr^leves et analyses par le laboratoire de Tusine. 

S . Les resultats de ces analyses ne seront done communiques a la fabrication 
qu'apres un nouveau temps mort, qui peut etre important (± 2 h). 
Modelisation du procede a Taide d'une fonction LAG 

Selon la methode de Tinvention, la modelisation dynamique d*un 
procede de synthese en continu est realisee en ayant recours h des hypo- 
10 theses de melanges parfaits et de delais purs. Les melanges parfaits sont 
mis en equation au moyen d'une fonction bien connue des ingenieurs, la 
fonction "LAG", ou filtre passe-bas (du ler ordre) ; cette fonction est 
lineaire et aisement programmable. 
Elle est definie comme suit : y = LAG (u , r) 
15 (dit "LAG de u durant t") comme etant la solution de Tequation 
differentielle 

dy 

dont les arguments u et r varient dans le temps. 

Cette equation peut se resoudre numeriquement (meme en temps 
reel) par une equation algebrique du ler ordre, ayant pour arguments les 
20 variables suivantes : 

- la periode d'echantillonnage "T" (ou temps ecoule depuis la demiere 
iteration) 

- le temps de sejour (ou "constante" de temps) "r" k Tinstant "t" 

- la grandeur d'etat "y" a I'instant precedent "t-T" 

25 - la grandeur de commande "u" k Tinstant present "t" 

(u et T represented en fait les vaieurs mesurees ou calculees a I'instant "t" 
sur les grandeurs "u" et "r", que Ton suppose avoir ete constantes dans tout 
I'intervalle precedent "T"), 

De preference, T est petit par rapport k t (par exemple au moins 10 
30 fois plus petit), de fagon a augmenter la precision du calcul. 

La solution de I'equation susmentionnee peut se calculer par exemple 
par la formule suivante : 

y(t) = y(t-T) . e + u(t) . ( I - e ""^^^(^^ ) 

ou encore, plus simplement (approximativement) : 
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y (t-T) + u(t) • ^ 



y(t) = 



1 + 



T(t) 



La modelisation du precede a Taide de la fonction LAG se base sur 
les theoremes suivants : 

Soit un reacteur parfaitement melange (CSTR), de volume Vj^. 
Divers constituants (reactifs ou inertes) ralimentent, panni lesquels le 
5 reactif "x" (d^bit entrant Fxj|^) ayant la propriety "Pxjjj" k rentnee. On 
mesure egaiement le debit de sortie "Fqxjt*' (soutirage). 
Thgorfeme 1 : Application de la m6thode LAG au calcul d*un bilan 
massique : 

A chaque instant, la masse "Mx^" d'un constituant "x" dans un 
10 reacteur parfaitement melange (CSTR) est egale au LAG du pnoduit du debit 
massique entrant "Fxjj^" par un temps "rj^", durant ce m6me temps "tx" : 
Mxj^ = LAG (Fxjj^ . Ty. , r^) (en kg) 

Le temps "r^** est le "temps de s6jour de x"; il vaut la masse du constituant 
Mx]^, divisee par la somme des debits massiques '*sortants" (quantite cons- 
15 omm6e par reaction "FRx", debit sortant du r&cteur "Fxquj", etc.). 

= Mxj^ / (Fxqut ^) 
Ce theoreme donne ainsi une methode exacte pour le calcul dynamique 
(meme en temps reel) des concentrations dans un reacteur parfaitement 
melange. En effet, la concentration "Cx^" du constituant "x" exprimee en 

20 kg/m^ vaut, si "Vj^" est le volume du reacteur ; 

Cxjj^ = Mxjj^ / Vjj^ (en kg/m^) 

Connaissant le debit volumique total "FVout" sortant du reacteur. 
on defmit le temps de sejour du reacteur "t^^" : 
TR = Vr / FVqut 

25 Dhs lors, le debit massique "Fxquy" du constituant "x" en sortie du 
reacteur (soutirage) vaut : 

On remarquera que si "x" est un inerte (ne subissant pas de reaction 
et ne sortant du reacteur que par soutirage), on a : 
30 rx = TR 

Par aiUeurs, dans le cas frequent oil la vitesse de reaction de "x" est 
proportionnelle a sa concentration Cxr, avec un facteur de proportionnalite 
R^, on a : FRx = R^ . Mxr et des lors 
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l/(Rx+ 1/tr) 

Thforfeme 2 : Application de la methode LAG au calcul d'une propriete 
d'un melange : 

Soit une propri6te "Px" d'un constltuant, repondant a la loi de 
5 melange lineaire : 

Pxi+2 = . Pxj + W2 , PX2 
oil wj et W2 sont les fractions massiques de propriete Pxj et Px2. (avec W| 
+ W2 = 1). 

A chaque instant la propriete Pxj^ dans un reacteur parfaitement 
10 melange (CSTR) est 6gale au LAG de la propriety k Tentrfc Pxq,^, durant 
un temps de sejour qui vaut le rapport de la masse M Xj^ du constituant dans 
le reacteur divise par le d€bit massique entrant (et/ou apparaissant par 
reaction) Fxjj^ : 

Pxr = LAG ( PxjN , Mxr / FxiN ) 
IS II est ainsi possible de tenir compte de la dynamique du proced^, et 

de recalculer en permanence ses constantes de temps. 

Comme expose prec6demment, les *'propriet6s" dont il est question 

ici peuvent dans certains cas avoir subi une transformation math^matique les 

rendant iindaires (par exemple, le logarithme de I'indice de fluidity en fondu 
20 d'un poiym^re peut etre consid6r6 comme r^ondant a une loi de melange 

lineaire). 

Principe de la regulation : 

Le module du proced6 de synthase etant etabli, il faut un algorittime 
qui calcule les param&tres n^cessaires a la regulation de ce proc&i^. 

25 La Figure 2 montre un sch6ma general simplifi6 du type de 

regulation avancee (Advanced Process Control ou "APC") adopte pour le 
proced6 de polymerisation d&rit plus haut. On voit que ce systeme de 
regulation comprend une cascade de deux algorithmes, cette cascade pilotant 
notamment les regulateurs PID des debits d'alimentation en reactifs. 

30 Les deux algorithmes en cascade, appeles algorithme maitre et algo- 

rithme esclave, sont tous deux des algorithmes dynamiques adaptatifs bases 
sur des modeles issus de la connaissance du processus (par opposition aux 
modeles empiriques), bases notamment sur les bilans de mati^re et la 
cinetique du procede regule. Us utilisent la fonction LAG pr6defmie. 

35 Sur la Figure 3, qui montre le principe du systeme de regulation 

dans le contexte du procede de polymerisation decrit plus haut, on voit 
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que : 

- ralgorithme maitre est base sur les equations caracteristiques des cata- 
lyseurs, c'est-k-dire des equations foumissant les proprietds du PE en 
fonction de la temperature de polymerisation et des concentrations des 

5 . reactifs dans le r6acteur; il foumit a ralgorithme esclave les consignes 
des concentrations en reactifs pour satisfaire les consignes pour les 
proprietes du PE; 

- ralgorithme esclave est base sur un bilan mati^re et la cinetique chimique 
des reactions; il foumit aux regulateurs des debits d'alimentation les 

10 consignes de dS)its de rdactifs ndcessaires pour satisfaire les consignes de 

concentrations imposees par ralgorithme maitre et la consigne d'allure de 
marche du procede. De preference, il calcule egalement un terme 
anticipatif ("feed-forward") pour la consigne de temperature, ameliorant 
la regulation de temf)erature iors des changements d'allure de marche. 
15 Ce type de regulation n'est d'une precision parfaite que si le modele 

est parfait et prend en compte toute les perturbations possibles. Cela n'est 
en general pas le cas. C*est pourquoi en general (voir Figure 4) on adapte 
en continu le module direct (et le modele inverse) en compaiant les 
predictions aux mesures des propri^t^s. Cette "adaptation" du modMe 
20 permet de maintenir sa precision en presence de perturbations non 
modelisees, et ainsi d'obtenir une regulation plus precise en toutes 
circonstances. 
L*algorithme esclave 

La Figure 6 iliustre le principe de Talgorithme esclave : 
25 1. un organe de prediction utilisant un module direct du procede predit, k 
paitir des mesures des debits alimentant le reacteur, les concentrations 
en reactifs et en polym&re; 

2. un organe d'adaptation compare les concentrations en ethylene (Et), 
hydrogene (Hy) et butene (Bt) mesuree par un analyseur (chroma- 

30 tographe a gaz) aux valeurs predites par le modele direct, de fagon a 

determiner trois parametres d' adaptation : 

- Tactivite specifique du catalyseur pour T&hylene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m*^ d'ethylene 

- Terreur de gain sur la mesure de d6bit d'hydrogene **KfHy" 
35 - la puret6 de Talimentation en butfene "KfBt" 

3. Torgane de contrdle calcule, a partir des consignes de concentration 
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calculees par ralgorithme maitre et de la consigne d'allure de marche, 
les consignes pour les debits d'alimentation des reactifs ; ces consignes 
de debit sont composees d'un feed-forward base sur le modele inverse, 
et d'un feed-back proportionnel k T&ait entre le modele direct et les 
5 consignes de concentration. 

Pour comprendre les calculs effectues par ralgorithme esclave, dans 
le cas d'un precede de synthese de polyethylene, U faut d'abord savoir qu'il 
est generalement admis que la vitesse de polymerisation "VitP^j" est 
proportionnelle : 

10 - k la concentration en Et non polymerise cEtjj (en kg/m^). 

- k la concentration en catalyseur actif dans le reacteur cCAjj^ (en kg/m^), 
et 

- a un facteur de proportionnalitd, fonction (difTicilement quantifiable) de 
la temperature, des concentrations en Hy, en Bt et en co-catalyseur, de la 

15 presence de poisons, etc. Ce terme est appele "activite catalytique" pour 

Tethylene Agj. Sauf dysfonctionnement majeur (poison, etc.), il varie 
raisonnablement peu en cours de campagne. 
VitPgt = Agj . cCAr . cEtg (kg/ra3;h) 

La quantite d'ethylene "FpEt" qui polymerise dans le reacteur a 
20 chaque unit6 de temps ("debit de polymerisation") vaut done, si Vj^ est le 
volume du reacteur : 

FpEt =VitPE^.Vj^ (kg/h) 

= Ag^ . cCAjj . MEtj^ 
ou MEtj^ est la masse d'ethylene en solution dans le reacteur (en kg) 
25 On sait par ailleurs que la vitesse d' incorporation de THy est environ 

100 fois plus lente que TR, et 10 fois plus lente pour le Bt. On en deduit : 
FpHy = Agj . cCAr . MHyj^ / 100 
FpBt = Ag^ . cCAj^ . MBtj^ / 10 

ou MHyj^ est la masse d'hydrogene en solution dans le reacteur (en kg) et 



30 MBtg^ est la masse de butene en solution dans le reacteur (en kg) 
Le mod^e esclave utilise maintenant les mesures suivantes : 

FEtjj,^ = ddbit d'alimentation en ethylene (monomere) (kg/h) 

FSvjj^= debit d'alimentation en solvant (hexane) (kg/h) 

FCAi|^= d6bit d 'alimentation en catalyseur (kg/h) 

35 FHyjj^= debit d'alimentation en hydrogfene (agent de transfert) (kg/h) 

FBtjj^ = debit d'alimentation en but&ne (co-monomere) (kg/h) 
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II utilise egalement les parametres d' adaptation suivants : 
^Et ^ "activite catalytique" pour I'ethylene 
KfHy = erreur de gain sur la mesure de Talimentation en hydrog^ne 
KfBt = puret6 de Talimentation en butane. 
5 Les calculs suivants sont effectues de facon sequentielle, dans I'ordre 

suivant, avec una frequence elevee (rintervalle de temps separant chaque 
iteration etant faible par rapport au plus petit des temps de sejour 7^). 

Le volume du reacteur etant consunt, ie debit volumique sortant 
equivaut au debit volumique entrant (fluides incompressibles). Le debit 
10 volumique soitant "FVqxjt ' P^"^ calcule comme la somme 

des d^its massiques entrants di vises par la masse volumique qu'ils ont dans 
le reacteur : 

FVoxjT = FSvji^/650 + FEtjj^/950 -f FBtjj^/eOO (m^/h) 
(ou les masses volumiques sont les suivantes : 650 kg/vn? pour le solvant, 

15 950 kglrn^ pour le polyethylene, 600 kg/m^ pour le butene). On fait ici 
rhypoth&se que tout Tethylene se transforme instantanement en poly- 
ethylene, et on neglige le debit d'hydrogene et de catalyseur (quelques kg). 

Le solvant est chimiquement ineite et ne sort du r&cteur que par 
soutirage du reacteur. On calcule sa masse "MSvj^" dans le reacteur en 

20 utilisant le theoreme 1 : 

Tj^ = Vj^ / FVOXJT (temps de sejour dans le reacteur) 

MSvj^ = LAG (FSvjj^ • ''R » ) O^g) 

Le catalyseur se desactive avec une constante de temps "kd" ; la 
masse "MCAj^" de catalyseur actif dans le reacteur est calculee comme 

25 suit : 

MCAr = LAG (FCAjjsj , tqj^ , tqj^ ) (kg) 
et done la concentration "cCAjj^" de catalyseur actif dans le reacteur vaut : 
cCAr = MCAr / Vj^ (kg/m^) 
30 L'6thylene "sort" du reacteur dans le debit de soutirage ainsi que 

dans la reaction de polymerisation. Sa masse dans le reacteur "MEtg^" est 
calcul6e comme suit : 

MEtR = LAG (FEtifj . , Tgt ) (kg) 
35 D'une fa^on similaire, on calcule la masse "brute" (non calibree) 

d'hydrogfcne dans le reacteur "MHyj^^" : 
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Tjjy = I / (1/tr + . cCAr / 100) (h) 

La masse "MHyj^", comg6e pour prendre en compte Teneur de 
gain sur la mesure de ralimentation en hydrog&ne, vaut : 
5 MHyR = BKHy . MHyj^^ (Kg) 

Le rapport Hy:Et dans le r^teur vaut done : 
HyEtjj = MHy^ / MEt^ 

D'une fa^on similaire, on calcule la masse de butene "brute" 

"MBtRAw"- 

10 ^Bt " ^ ^ (^^^R ^Et • ^^^R ' 

MBtj^W = (™^IN • ^Bt • ^Bt ) 

La masse "MBtj^", corrigee pour prendre en compte la purete de 

ralimentation en butane, vaut : 

MBtR = KfBt . MBtj^w (kg) 
15 Le rapport Bt:Et dans le reacteur vaut done : 

BtEtj^ = MBtj^ / MEtj^ 

On a vu que le d6bit de polymerisation "FpEt" (allure de marche 
instantan6e) vaut : 

FpEt = Agt . cCAr . MEtjj (kg/h) 
20 Le polymfere etant inerte et ne subissant pas de reaction, sa masse 

dans le reacteur "MPEp^" vaut : 

MPEjj = LAG ( FpEt . tr, tr) (kg) 
Le debit de polymere sortant du reacteur "FPEqut" ^aut des lors : 
FPEqut = MPEr / (kg/h) 
25 Adaptation d u modele esclave : 

Le bloc d'adaptation utilise un analyseur (par exemple un chroma- 
tographe k gaz "GC") pour obtenir les mesures des concentrations dans le 
reacteur en ethylene "cRq^", en hydrogene "cHyQ^" en butene 
"cBtQQ" (exprimees par exemple en kg/m^). Ces mesures sont comparees 
30 aux valeurs predites par le modfele direct, afm de determiner les trois 
param^tres d'adaptation suivants : 

- ractivit6 sp&ifique du catalyseur pour I'^thylene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m^ d'dthylfene 

- Terreur de gain sur la mesure de debit d'hydrogene "KfHy" 
35 - la purete de ralimentation en butene "KfBt" 

Le chromatogiaphe a gaz foumit des mesures echantillonnees avec 
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un retard d'environ 6 miiiutes. 

Le calcul de I'activite specifique pour I'ethylene "AEt" decoule des 
6galites suivantes : 

cBLqq = cEtj^ (kg/m^) 
5 = MEtR / Vj^ 

= LAG (FEtjN . , ) / Vr 
= LAG (FEtjN / (I/XR + Agt . cCAr) , ) / 
Agj" = ( FEtnq / LEAD ( cRqC ' . ret ) - 1/tr ) / cCAr 

La mesure cELqq ^tant 6chantillonnee et bruits, on ne souhaite pas 
10 la trouver dans un LEAD; dhs lors, on utilise plutdt : 
A^' = LAG (Aej" , 

« LAG (FEtjj^/cCAR, tj^ I {cEXqqM^ - 1 / LAG (tr.cCAr, Tj^) 
On tient compte du retard de 6 minutes sur la mesure de cEtQ^;; en 
introduisant deux LAGs en serie de 3 minutes chacun sur les valeurs du 
15 modele, et on obtient la formule finale : 

A^t = LAG ( FEtjj^ / cCAr , , 3/60 , 3/60 ) / (cEtQ^ . Vr) 
- 1 / LAG ( TR . cCAr , rgt , 3/60 . 3/60 ) 
Le calcul du gain du debit d'hydrogene "KfHy** d&ouie des ^galitds 
suivantes : 

20 cHyQc = cHyR (kg/m^) 
= MHyR / Vr 
= KfHy . MHyRAW ' ^r 
KfHy= cHyQc • Vr / MHyRAW 

On tient compte du retard de 6 minutes sur la mesure de cHyQ(^ en 
25 introduisant deux LAGs en serie de 3 minutes chacun sur la valeur du 
modMe, et on obtient la formule finale : 

KfHy = cHyQc . Vr / LAG ( MHyR^w » 3/60, 3/60 ) 
D'une fa^on similaiie a KfHy, on calcule le parametre de correction 
pour la purete du butene "KfBt" : 
30 KlHy = cBtQc . Vr / LAG ( MBtR^w > 3/60, 3/60 ) 
L'algorithme de regulation a pour entrees : 
- les consignes de concentration calculees par Talgorithme maitre; plus 
pr6cisement les consignes pour les rapports de concentration cHyR/cRR 
"HyEtsp" et cBtR/cEtR "BtRgp" (en kg/kg) 
35 - la consigne pour Tallure de marche du procede FpEtgp, fixee par 
Toperateur 
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- la consigne pour la concentration en ethylene cEtgp, fix6e par 
rop6rateur, 

- les concentrations calcul6es par le module. 

II calcule les consignes pour les debits d'alimentation des reactifs 
5 FEtgp, FCAgp, FHygp et FBtgp. Differents algorithmes peuvent etre 
utilises, parmi lesquels le MBPC (Model Based Predictive Control). lis 
peuvent en general se decomposer en un feed-forward base sur le module 
inverse, et un feed-back proportionnel a Tecart entre le modele direct et les 
consignes de concentration. 
10 Regulation de ralimentation en Ethylene : 

Feed-forward : valeur pour maintenir la concentration actuelle, bas6e sur 
rinversion de la valeur stationnaire de Tequation suivante : 

MEtR = LAG (FEtiN • ^Et ' ^Et > " (^S) 

FEtpF = MEtR / 

15 Feed-back : proportionnel a Tecart entre la consigne cEtgp, et le modele 
FEtpB = 5 . (cRgp . Vr - MEtR) 
Consigne : 

FEtgp = FEtpp + FEtpB 
R6pulation de ralimentation en catalyseur : 
20 Feed-forward : valeur pour maintenir la concentration actuelle, basee 

sur rinversion de la valeur stationnaire de Tdquation suivante : 

MCAr = LAG (FCAiN • ^CA ^ ^CA ) (^g) 
FCApp = MCAR/rcA 
Feed-back : proportionnel a Tecart entre la consigne FpEtgp, et le modele, 
25 selon la formule suivante : 

FpR = . MEtR , MCAr / Vr 
FCApB = 5 . (FpEtgp / (Agj . MEtj^ / Vj^) - MCAr) 
Consigne : 

FCAgp = FCApp + FCApB 
30 Regulation de ralimentation en hvdropene : 

Feed-forward : valeur pour maintenir la concentration actuelle, basee 
sur rinversion de la valeur stationnaire de I'dquation suivante : 

MHyR = LAG (FHyi^ . r^y , ryy ) (kg) 

FHypF = MHyjjAw / ^Hy 
35 Feed-back : proportionnel a recart entre la consigne de rapport HyEtgp, et 
le modele. 



I 
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FHypB = 5 . (HyEtsp • MEt^ - MHy^) 
Consigne : 

FHysp = FHypp + FHypg 
Repulation de ralimentation en butfene (similaire a Thydrogene) : 
5 . FBtpp = MBtRAW ' ^Bt 

FBtpB = 5 . (BtEtsp . MEtj^ - MBtj^) 

FBtgp = FBtpp + FBtpg 

Les equations qui precedent resument les equations de I'algorithme 
esclave. Elles sont executees toutes les 10 secondes par le systeme nume- 
10 rique de controle et de commande (SNCC). 
L'algorithme maitre 

La Figure 5 illustre le principe de I'algorithme maitre : 

1. son organe de prediction (base sur un modele direct) predit les 
proprietes principales du polymere (MI et MVS); il utilise pour cela la 

15 mesure de la temperature de polymerisation, les predictions pour les 

concentrations dans le reacteur, foumies par le module esclave, et les 
temps de sdjour du PE dans les diffdrents appareils ; 

2. son organe d'adaptation compare les mesures de MI et de MVS 
effectu^ (a la sortie du secheur egalement), soit par le laboratoire de 

20 mesure toute les 2 h, soit par un analyseur en continu, aux valeurs 

predites par le modele direct, de fa?on a determiner les 2 parametres 
d'adaptation, qui sont des parametres correctifs, multiplicatif pour le MI 
et additif pour la MVS ; 

3. son organe de controle (base sur un module inverse) calcule, a partir des 
25 consignes de MI et de MVS foumie par Toperateur, les consignes pour 

les concentrations dans le reacteur (n^orts Hy:Et et Bt:Et). Comme 
pour Talgorithme esclave, ce calcul est compose d'un feed-forward base 
sur le modele direct, et d'un feed-back proportionnel a Tecart entre le 
modele direct et les consignes de Toperateur. 
30 Pour un catalyseur donne, les proprietes des resines en regime stable 

sont des fonctions de la temperature de polymerisation et des concentrations 
des reactifs. Parmi les differentes Equations statiques renseignees dans la 
litt6rature, on a choisi les equations suivantes : 
log(MI) = ao + aj . -f a2 . log (Hy/Et) + a3 . Bt/Et 
35 MVS = bg + bj . T** + b2 . (Bt/Et)'^^ + b4 . log(MI) 

Les parametres ag a a^ et bQ a b^ sont obtenus par identification en 
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regime stable, pour plusieurs resines fabriquees avec le meme catalyseur. 

Par aiUeurs, les differents appareils rencontres par le polyethylene 
jusqu'au moment ou ses proprietes sont mesurees (neacteur, stripper, 
centrifiigeuse puis secheur), peuvent tous etre assimiles en premiere appro- 
5 ximation a des reacteurs parfaitement melanges. 

L'algorithme maitre utilise comme entrees les mesures suivantes : 
Tt^ = temperature au rdacteur ("C) 
Vstp = volume liquide dans le stripper (obtenu par mesure de niveau) (m*') 
^MES ^ n^esure de MI (melt-index) 
10 ^^Sj^^jgg = mesure de MVS (Masse Volumique Specifique) 
ainsi que les calculs suivants realises par Talgorithme esclave : 
FpEt = production instantanee de polymere (allure de marche) (kg/h) 



FPEqutj, = debit de PE sortant du reacteur ' (kg/h) 
MPEg^ = masse de PE dans le reacteur (kg) 

15 HyEtj^ = rapport Hy sur Et dans le r^cteur (kg/kg) 

BtEtp^ == rapport Bt sur Et dans le r^cteur (kg/kg) 



Les valeurs instantanees brutes (avant adaptation) de la MVS et du 
logarithme du MI ("IMI") sont calculees par : 
MVSiNs = bo + bi . Tr + b2 . (BtEtR)b3 + b4 . IMIijjs 
20 JMIijvjs = aQ + aj . Tj^ + a2 . log (HyRj^) + a3 . BtEtj^ 

En utilisant le theoreme 2, on calcule les proprietes moyennes brutes a 
la sortie du reacteur : 
IMIr = LAG (IMIjj^g, MPEj^ / FpEt) 
MVSr = LAG (MVSn^g, MPEr / FpEt) 
25 En effet : 

- les proprietes IMI et MVS r^ndent assez bien a une ioi de melanges 
lineaire 

« le reacteur boucle peut etre assimile a un reacteur parfaitement melange 

- le d6bit massique de PE "entrant" (apparaissant) dans le reacteur est 
30 effectivement FpEt, la quantite de PE qui polymerise a chaque instant 

(allure de marche). 

Proprietes brutes a la mesure : Sachant qu'on a dans le stripper 
environ 500 kg de PE par m^, et si Ton prend Thypothese que le stripper 
est un reacteur parfaitement melange, on calcule les proprietes brutes en 
35 sortie du stripper comme suit : 

IMIstp = LAG (IMIr, 500 , Vstp / FPEqut) 
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MVSstp = LAG (MVSr, 500 . Vstp / FPEqut) 

Le temps de sejour dans la centrifugeuse 6tant tres court, on pent 
rignorer. 

Le secheur est a lit fluide; il contient en permanence environ 
5 1. 400 kg de PE. On peut faire I'hypothese que le niveau dans le stripper 

change peu, et que le debit qui en sort egale celui qui y rentre. Des lors, le 

debit de PE entrant dans le secheur vaut FPEoxjt- On a alors en sortie du 

secheur, k Tendroit oil T&hantillon est preleve pour la mesure de propriety, 

les valeurs brutes suivantes : 
10 IMIsh = LAG (IMIstp, 1400 / FPEqxjt) 

MVSsh = LAG (MVSstp. 1400 / FPEqut) 

Les proprietes apr&s adaptation sont obtenues en faisant intervenir les 

parametres d'adaptation kMI (parametre multiplicatif) et kMVS (parametre 

additif); les proprietes aprfes adaptation sortie r^cteur, stripper et secheur 
IS valent done : 

Mlr^ = kMI • lO'^^ 

MVSrc = kMVS + MVSr 

MIstpc = kMI . lO^^stp 

MVSstpc = KMVS + MVSstp 
20 MIshc = kMI . lO^^sh 

MVSshc = kMVS + MVSsh 

Adaptation de Talgorithme maitre : 

Les mesures de proprietes prennent un certain temps a etre realisees 

(± 5 min. si analyseur en ligne, ± 1 h si effectuees par le laboratoire). 
2S Pour pouvoir calculer les param6tres d'ad^q^tation, il faut done 

resynchroniser (recaler dans le temps) les predictions brutes du module avec 

les mesures. Cela peut par exemple se faire au moyen d'un registre a 

decalage (appele ici "fonction DELAY") : 

IMIj^EL E>ELAY (IMIsh, rj^) 
30 MVSi^gL = DELAY (MVSsh, Tj^ys) 

avec : rj^ et tj^vS ~ i 5 min. ou ± 1 h, suivant que la mesure est 

realisee par un analyseur continu ou par le laboratoire. 

A cbaque nouvelle mesure de MI ou de MVS, on recalcule le 

parametre d'adaptation brut kMI* ou kMVS\ en comparant la valeur 
35 modele brute resynchronisee a la valeur mesuree : 

kMI' = log (M1meS> ' "^DEL 
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kMVS' = MVSmes - MVSdel 

Ces valeurs brutes sont filtr&s afin d*att^nuer les reactions 
precipitees que les peiturbations (bruit) eventuelles des mesures infligeraient 
au procede : 
5 kMI = LAG ( kMI', ± 1 h) 

kMVS = LAG ( kMVS\ ± 1 h) • 
Orpane de controle : 

L'organe de controle a pour consignes les valeur MIgp et MVSgp, 
entrees par Toperateur. II calcule les consignes pour les rapports des 
10 concentrations dans le reacteur HyEt^^p et BtEt^p necessaire pour obtenir 
rapidement les proprietes voulues MI^p et MVS^p. Ce calcul se fait en 2 
etapes : 

1. Torgane de controle calcule, a partir des consignes MIgp et MVSgp 
foumies par I'opdrateur, et des valeurs aprfes adaptation du MI et de la 

15 MVS dans les differents appareils, les consignes Mligp et MVSigp pour 

la production instantanee. Ces consignes instantanees sont composees 
d'un feed-forward, et d'un feed-back propoitionnel k I'&art entre le 
modMe direct et les consignes de Topdrateun 

2. les consignes pour les rapports de concentration HyEt^p et BtEtj^p sont 
20 alors calculees en inversant Tequation statique utilisee plus haut pour le 

calcul de la valeur instantanee du MI et de la MVS. 

Consignes pour les proprietes instantanees : On compare les 
proprietes sortie secheur aux consignes de proprietes, pour determiner les 
consignes voulues pour les propriet6s a la sortie du stripper (on neglige la 
25 centriAigeuse) : 

MIstpsp = 10 ( log(MISP) + 0.1 . (log(MISP) - log(MIshC))) 

MVSstpsp = MVSgp + 0.1 . (MVSgp - MYSsh^) 

De mfeme, a partir de Tecart entre ces consignes sortie stripper et 

des valeurs calibrees au stripper, on calcule les consignes voulues en sortie 

30 du reacteur : 

Mir«^ r= 10 ( Iog(MIstpSP) + 0.5 . (log(MIstpSP) - log(MIstpC))) 

MVSrgp = MVSstpsp + 0,5 . (MVSstpgp - MVSstp^) 

Enfin, k partir de Tecart entre ces consignes sortie reacteur et des 
valeurs calibrees correspondantes, on calcule les consignes voulues pour la 
35 production instantanee : 

MlicD = 10 ( log(MIrSP) -f 2 . (log(MIrSP) - log(MIrC))) 
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MVSigp = MVSrgp + 2 . (MVSrsp " MVSfq) 

Consignes pour les rapports de concentrations : Les consignes pour 
les rapports de concentration HyEt^p et BtEt<jp s'obtiennent en inversant 
r equation statique utilisee plus haut pour le calcul de la valeur instantanee 
5 du MI et de la MVS, en substituant aux termes MI et MVS les consignes 
voulues pour la production instantanee et en appliquant le parametre 
d'adaptation. 

Partant de : 

log (Mligp / kMI) = aQ + aj . Tj^ + a2 . log (HyEtgp) -h sl^ - BtEtj^ 
10 MVSigp - kMVS = bg + bj . Tr + b2 . (BtEtsp)''^ + b4 . IMIn^g 
on obtient : 

a2 . log (HyEtgp) = log (Mligp / kMI) - (Hq + aj . Tjj + 33 . BtEtj^) 
ce qui donne : 

HyEtgp = 10 ( ( (MliSP / kMI) - aO - al , TR - a3 . BtEtR ) / a2 ) 
15 et : 

b2 . (BtEtsp)*^^ = MVSigp . kMVS - ( bg + bj . + b4 . IMIjj^s ) 
ce qui donne : 

BtEtgp = ( ( MVSigp - kMVS - bg - . Tr - b4 , IMIn^g ) / b2 ) 

Les Equations qui precedent resument les equations de Talgorithnie 
20 maitre. EUes sont executees toutes les 30 secondes par le SNCC. 

Avec ce proced6, il est possible de piloter la polymerisation avec une 

grande precision. En particulier : 

- les proprietes controlees (MI et MVS) sont maintenues au plus pres des 
valeurs voulues, avec une dispersion minime 

25 - les changements de qualite (et done des proprietes MI et MVS) sont 
realises avec rapidite et precision 
" les demarrages et arrets de la polymerisation, ainsi que les changements 
d 'allure de marche du procede, sont realises d'une fa^on acceleree, tout 
en maintenant le MI et la MVS fort proches des valeurs voulues. 
30 Bien que la methode de regulation selon T invention ait ete presentee 

k Taide d'un proced^ de synthese de polyethylene par polymerisation en 
continu d*ethylene, il est entendu que cette methode de regulation sera de 
faQon g6n6rale efflcace pour d^autres precedes de synthese, et en particulier 
pour les procedes presentant Tune ou plusieurs des caracteristiques 
35 suivantes : 

- une regulation muhivariable est necessaire parce que plusieurs variables 
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influent sur Tensemble des proprietes a regler; 

- la dynamique du processus est lente : melanges en serie, temps morts 
■importants; 

- les mesures des proprietes sont &hantillonnees avec une faible frequence, 
5 et/ou bruitees; 

- la regulation doit etre dynamique, c'est-a-dire valable quelle que soit 
Tallure de marche du procede, ainsi que durant les transitions d'allure de 
marche et de qualite (proprietes) du produit k synthetiser; 

- 11 est interessant d'estimer certaines variables non directement mesurees. 
10 Pour pouvoir etre facilement mise en oeuvre avec les techniques 

presentees, il suffit que : 

- les equations statiques du processus soient connues (souvent, elles le sont 
au moins dans une certaine mesurc, sinon le procedd ne pourrait etre 
maitrise); 

15 - la dynamique du processus puisse etre approchfe par des melanges 
parfaits et des temps morts; 

- les mesures necessaires soient disponibles et d'une quality suffisante (en 
particulier, les debits des reactifs et les debits traversant les reservoirs 
concern^). 

20 L'utilisation particuliere de la fonction LAG decrite ci-dessus, 

notamment dans les theoremes 1 et 2, pent naturellement etre etendue a des 
methodes de regulation basees sur une structure autre que celle qui a 6i6 
exposee, qui comprend un algorithme maitre et un algorithme esclave 
distincts. EUe peut par exemple fetre appliquee dans une methode de 

25 regulation ne comportant qu'un algorithme unique. 
Exemple? : 

8 essais de synthese de polyethylene (PE) de 4 types differents 
(defmis par leur MI, MVS, ...) ont ete realises, respectivement en utilisant 
une methode de regulation classique et en utilisant la methode de 
30 rinvention. Le tableau ci-dessous resume les constatations qui ont ete faites 
sur la base de nombrcuses mesures de I'indice de fluidite en fondu (meh 
index) des 8 polymeres obtehus. Cpk d&igne Tindice de capabilite centre 
du procede. 
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Type 
. dePE 


Regulation classique 


Regulation selon I'invention 


Ecart-type 


Cpk 


Ecart-type 


Cpk 


1 


0.127 


0.909 


0.059 


2.202 


2 




0.61 




2.0 


3 




0.48 




1.88 


4 




0.64 




1.09 



On constate que Tindice de capabilite Cpk est plus que double grace 
a Tutilisation de la methode de Tinvention, ce qui indique que les proprietes 
sont environ deux fois moins dispersees et/ou mieux centrees par rapport 
aux consignes imposees. 
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Tableau des abreviations utilisees ' 

Hj parametres pour Tequation statique du MI (I = 0 a 3) 

bj parametres pour I'equation statique de la MVS (I = 0 a 4) 

Agj activite catalytique pour Tethylfene (m^.kg'^h"^) 

5 cxqq concentration de "x" obtenue par la mesure de Panalyseur (kg/m^) 

cxjj concentration de "x" dans le reacteur (kg/m^) 

cxgp consigne pour la concentration de "x" dans le i^teur (kg/m^) 

Fpx debit massique de polymerisation de "x" (allure de marche) (kg/h) 

FVquj debit volumique sortani du reacteur (tn^/h) 

10 Fxjj^ debit massique de "x" entrant (kg/h) 

Fxquj debit massique de "x" sortant (kg/h) 

kd constante de desactivation du catalyseur (i/h) 

KfBi parametre correctif (d'adaptation) pour le butbne 

KfHy panun&tre correctif (d'adaptation) pour Fhydrogene 

15 kMI parametre correctif (d'adaptation) pour le MI 

kMVS parametre correctif (d'adaptation) pour la MVS 

LAG ( , ) fonction filtre passe-bas du ler ordre 

^MES mesure du MI (Melt Index) (indice de fluidite en fondu) 

MIy MI (Melt Index) brut (non calibr6) dans "y" 

20 Mly^ MI calibre (avec adaptation) dans "y" 

Mlygp consigne de MI (calibree) pour "y" 

MVS|^£5 mesure de la MVS (Masse Volumique Standard) 

MVSy MVS (Masse Volumique Standard) brute (non calibre) dans "y" 

MVSy^ MVS calibree (avec adaptation) dans "y" 

25 MVSygp consigne de MVS (calibree) pour "y" 

^^RAW masse de "x" brute (non calibre) dans le reacteur (kg) 

Mxy masse de "x" calibree (avec adaptation) dans "y" (kg) 
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Rj^ reactivite de X dans le reacteur 

Vy Volume de "y" (m^) 

temps de sejour dans le reacteur (h) 
r^' temps de sejour pour "x" dans le reacteur (h) 

5 "x" peut representer les constituants suivants : 
Bt Butene 
CA Catalyseur 
Et Ethylene 
Hy Hydrogene 
10 Sv Solvant 

"y" peut representer les appareils suivants : 
r Reacteur de polymerisation 

stp Stripper 
sh Secheur (a lit fluide) 
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Le^ende des figures 

10 reacteur de polymerisation 

11 alimentation en rdactifs (matiercs premieres), caialyseur, solvant 

12 circuit de refroidissement 
5 14 conduite de soutirage 

16 stripper 

18 condenseur 

20 chromatographe en phase gazeuse 

22 centrifugeuse 

10 24 secheur a Ut fluide 

25 solvant et reactifs a recycler 

26 polyethylene 

27 reactifs a recycler 

28 solvant k recycler 

15 30 consignes de propriete du polymere 

31 consigne d'allure de marche du procede 

32 algorithme maitre 

33 algorithme esclave 

34 consignes de concentration 
20 35 consignes de debits entrants 

36 regulation (PID) des debits 

37 mesures 

38 regulation de temperature 

39 feed-forward de temperature 

25 40 dynamique de la polymerisation : cinetique chimique et bilan matiere 

41 simulations de Failure de marche et des rapports Hy/Et et Bt/Et 

42 mesures de temperature, de debits et de concentrations 

43 procede regule 

44 analyse d'un echantillon du produit synthetise par le procede 
30 51 mesure de grandeurs liees au deroulement du procede 

52 mesure de proprietes du polymere 

53 modele direct ; prediction des proprietes i la mesure 

54 comparaison : correction du modele (adaptation) 

55 algorithme de regulation, base sur le modele inverse (feed-forward + 
35 . feed-back) 
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56 consignes pour les grandeurs liees au procede 

57 mesures el consignes de debits entrants . 

58 mesure de la temperature et prediction des concentrations 

59 modele direct : equations des proprietes en fonction des concentrtations 
5 60 prediction des proprietes 

61 consignes pour les concentrations dans le reacteur 

62 mesure des debits entrants 

63 mesure des concentrations dans le reacteur 

64 modele direct ; pr&iiction des concentrations sur base du bilan matiere 
10 65 pr^ictions de concentration 

66 comparaison ; calcul des panunetres d'adaptation 



PCT/EP96/02408 



-43-. 

RE VENDIC ATIQNS 

1 - Methode de regulation d'un proc6d6 de synthase d'au moins un 
produit chimique dans un appareillage comprenant au moins un reacteur (R) 
pouvant etre assimile h un reacteur parfaitement melange, dans laquelle une 
ou plusieurs grandeurs de commande (GC) permettent d'agir sur le 
d6roulement du precede en vue qu'une ou plusieurs grandeurs liees aux 
propriet^s du produit et/ou au deroulement du procede, appelees grandeurs 
regi&s (GR), soient 6gales a des consignes corespondantes (Cq|^), ladite 
methode comprenant les etapes suivantes : 

(a) entree de consignes concemant les grandeurs reglees (Cq^) ; 

(b) calcul, au moyen d'un organe de prediction (OP), de predictions des 
grandeurs reglees (Pqr), sur la base de mesures des grandeurs de 
commande du procede (Mq^) ; 

(c) utilisation d'un organe de controle (OC) pour calculer des consignes 
des grandeurs de commande du proc^e (Cqq), sur base des 
consignes (Cqj^) et des pr&Iictions (Pqr) des grandeurs reglees ; 

(d) transmission des consignes des grandeurs de commande du procede 
(Cqq) k des actionneurs, ou a des organes de reglage controlant des 
actionneurs, afm d'agir sur le deroulement du procede ; 

dans laquelle T organe de prediction (OP) est bas6 sur un modele 
mathematique du procdde, appele modele direct (M), et est con?u de telle 
manifere qu'on predit la masse Mjq^ d*au moins un constituant (X) dans le 
reacteur (R) par Tdquation : 

MxR = LAG (FxRin • » > 

dans laquelle : 

■ ^XRin ^^^^^ massique du constituant X entrant dans le reacteur R ; 

- Tj^ est le temps de sejour de x dans le reacteur, qui vaut - 

^X = / (E Fxdis) 

dans lequel : 

^XR d^signe la demiere valeur calculee de la masse du constituant 
X presente dans le reacteur R ; 

E Fxdis designe la somme de tous les debits massiques Fxdis avec 
lesquels le constituant X dispardt du r^cteur R, notamment par 
reaction et/ou par sortie du reacteur ; 

- la fonction y = LAG (u , r) est la solution de Inequation differentielle 



I 
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u = T ' 4^ + y 
dt ^ 



calculee avec la valeur instantanee de u et de r, ainsi qu'avec la demiferc 
valeur de y calculi. 

2 - Methode de regulation selon la revendication 1 , dans laquelle la 
consigne d'au moins une grandeur reglee (Cqjj) est corrigee sur la base de 
I'ecart entre la mesure (Mqj^) et la prediction (Pqr) de cette grandeur 
r6glee, de fagon k ce que la regulation soit efficace meme en presence d'une 
erreur dans la prediction de cette grandeur reglfe (Pgr). 

3 - Methode de regulation selon la revendication 1, dans laquelle le 
modele (M) du proced6 est periodiquement adapte sur la base de Tecart 
entre les predictions (Pqr) et les mesures (Mqr) des grandeurs reglees, de 
fagon a ce que le modele du procede foumisse des predictions des grandeurs 
reglees (Pqr) aussi proches que possible des mesures de ces grandeurs 
(Mqr). 

4 - Methode de regulation selon la revendication 3, dans laquelle les 
IS mesures (Mqr) des grandeurs r^gldes interviennent uniquement dans 

Tadaptation eventuelle du modele du procede, et n'interviennent pas 
directement dans le calcul des consignes des grandeurs de commande du 
procede (Cqq). 

5 - Methode de regulation selon Tune des revendications 114, 
20 appliquee k un procddd de polymerisation, comprenant Tune ou plusieurs 

des 6tapes supplementaires suivantes : 

- calcul d*une consigne de temperature dans le reacteur en fonction d*une 
ou plusieurs consignes de proprietes du produit ; et transmission de cette 
consigne de temperature a un ou plusieurs actionneurs permettant de 

25 modifler la temperature dans le reacteur ; 

- calcul d'un bilan thermique pour le reacteur, sur base notamment de 
mesures de temperature; utilisation de ce bilan thermique de fagon a 
determiner la quantite de polymere synthetise par unite de temps et/ou la 
productivite du catalyseur et/ou la concentration d'au moins un reactif 

30 dans le reacteur ; 

- calcul de la quantite de chaleur produite par la polymerisation, grace a un 
calcul de la quantite du ou des reactifs qui polymerisent ; determination 
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par ce biais de la quantite la chaleur qu'il faut ajouter ou evacuer pour 
maintenir la temperature du reacteur ; utilisation du resultat dudit calcul 
pour ameiiorer la regulation de temperature, de fa^on a respecter au 
raieux la consigne de temperature, notamment en cas de modifications de 
S . Failure de marche. 

6 - Mdthode de regulation selon Tune des revendications 
pr&edentes, dans laquelle la propriete Pxjj d'un constituant "x" dans le 
reacteur R, assimile a un reaaeur parfaitement melange, est calculfe comme 
suit : 

10 Pxr = LAG (PxjN , Mxr / FxjN ) 

ou "Px" est une propriete d'un constituant "x", repondant de maniere 
substantielle a la loi de melange lineaire Pxj^2 Wj . Pxj + W2 . Px2 , 

et W2 etant les proportions massiques de deux fractions 1 et 2 de 
propriete Pxj et Px2 que Ton melange ; 

15 Pxj ^ propriety de x k sa sortie du r&cteur, aprfes melange ; 

Pxjrj^ est la propriete du constituant "x" a son entree dans le r&cteur R ; 

Mx^ est la masse du constituant x dans le reacteur R ; 

Fxq^ est le ddbit massique du constituant x entrant dans le r^cteur R. 

7 - Methode de regulation selon Tune des revendications precedentes, 
20 comprenant les etapes suivantes : 

- entree de consignes relatives a une ou plusieurs proprietds du produit a 
synthetiser, dans un algorithme maitre ; 

- entree de la consigne d*allure de marche du procede, dans un algorithme 
esclave; 

25 - calcul de consignes de concentration des constituants dans le reacteur a 

Taide de Talgorithme maitre, en fonction notamment des consignes et des 
mesures de propriety du produit ainsi que de mesures.ou de predictions 
des concentrations des diff6rents constituants dans le r&icteur ; 

- transmission des consignes de concentration calculees par Talgorithme 
30 maitre comme grandeurs d'entree a I'algorithme esclave; 

- calcul des consignes de debit des constituants entrant dans le reacteur, k 
Taide de Talgorithme esclave, en fonction notamment de la consigne 
d'allure de marche du procede, des consignes de concentration, et de 
mesures de debit des constituants entrant dans le reacteur, et 

35 - transmission des consignes de debit calculees k I 'aide de Talgorithme 
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esclave a un ou plusieurs actionneurs afin de regler les debits des 

constituants entrant dans le leacteur, 
dans laquelle ralgorithme maitre et/ou l*algorithme esclave sont mis en 
oeuvre confonnenient k Tune des revendications precedentes. 

8 - Methode de regulation selon la revendication 7, caiacterisee en 
ce que ralgorithme maitre comprend : 

- un organe de prediction, base sur un modele direct du precede permettant 
de foumir une prediction des proprietes du produii synthetise en fonction 
de mesures et/ou de predictions des concentrations des constituants ; 

- un oi;gane d'adaptation comparant les predictions de proprietds calculdes 
par Torgane de prediction a des valeurs effectivement mesurees sur le 
produit synthetise et deduisant de cette comparaison des parametres 
d'adaptation, iesdits parametres d'adaptation intervenant comme 
grandeurs d'entrte supplementaires dans ledit organe de prediction de 
ralgorithme maitre ; et 

- un oigane de contrdle, bas6 sur un module inverse du procede, pour 
calculer, en fonction des consignes et des predictions de propri6t6s du 
produit k synthetiser, des consignes de concentrations pour ralgorithme 
esclave, Iesdits parametres d'adaptation intervenant egalement comme 
grandeurs d'entree supplementaires dans ledit organe de controle. 

9 - Mdthode de regulation selon Tune quelconque des revendications 
7 ou 8, dans laquelle ralgorithme esclave comprend : 

- un organe de prediction, base sur un modele direct du procede permettant 
de foumir une prediction des concentrations d'un ou plusieurs des 
constituants sur base d'un biian de mati^re dans le reacteur ; 

- un organe d'adaptation, comparant les predictions de concentrations 
caiculdes par le modele direct k des mesures de concentration et 
deduisant de cette comparaison des parametres d'adaptation, Iesdits 
parametres d'adaptation intervenant comme grandeurs d'entree 
supplementaires dans ledit organe de prediction de ralgorithme esclave ; 
et 

- un organe de controle, base sur un modfele inverse du procede, pour 
calculer, en fonction de la consigne d'allure de marche, des consignes de 
concentration calcuiees par rorgane de contrdle de I'algorithme maitre et 
des predictions de concentration calcuiees par rorgane de prediction de 
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ralgorithme esclave, les consignes pour les debits entrant dans le 
reacteur, lesdits param&tres d'adaptation intervenant comme grandeurs 
d'entree supplementaires dans ledit organe de contrdle de I'algorithnie 
esclave. 

5 • 10 - Methode de regulation selon Tune des revendications 

pr&6dentes, appliquee a la regulation de la synthese en continu de 
polyethylene par polymerisation d'ethylene dans au moins un reacteur, les 
reactifs comprenant Tdthylene, Thydrogene et/ou un comonomere optionnel, 
la reaction de polymerisation ayant lieu en presence d*un catalyseur et une 
10 partie du contenu du reacteur etant preleve en permanence ou par 
intermittence. 

11 - Methode de regulation selon les revendications 9 et 10, dans 
laqueUe Toigane d'adaptation de Talgorithme esclave compare les mesures 
des concentrations en ethylene (Et), hydrogfene (Hy) et/ou comonomfere 

15 optionnel (Bt) aux valeurs predites par Toigane de prediction de 

Talgorithme esclave, de fagon k determiner au moins Tun des paramfetres 
d'adaptation suivants : 

a) I'activitd specifique du catalyseur pour Tethylfene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m^ d'ethylene ; 
20 b) I'erreur de gain sur la mesure de debit d'hydrogene "KfHy" ; 
c) la puret6 de Talimentation en comonomere "KfBt". 

12 - Methode de regulation selon Tune des revendications 1 k 9, 
appliquee a la regulation de la synthase en continu de polypropylene par 
polymerisation de propylene dans au moins un reacteur, les reactifs 

25 comprenant le propylene, Thydrogfene et/ou un comonomfere optionnel, la 
reaction de polymerisation ayant lieu en presence d'un catalyseur et une 
paitie du contenu du reacteur etant preleve en permanence ou par 
intermittence, 

13 - Methode de regulation selon les revendications 9 et 12, dans 
30 laquelle Torgane d'adaptation de Talgorithme esclave compare les mesures 

des concentrations en propylene (Pe), hydrogene (Hy) et/ou comonomere 
optionnel (Et) aux valeurs predites par Torgane de prediction de 
Talgorithme esclave, de fagon a determiner au moins Tun des parametres 
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d'adq^tation suivants : 

a) Tactivite specifique du catalyseur pour le propylene "APe", en kg/h de 
polypropylene par kg de catalyseur et par kg/m-^ de propylene ; 

b) rerreur de gain sur la mesure de debit d'hydrogfene "KfHy" ; 
5 c) la piirete de T alimentation en comonomfere "KfEt". 

14 - M6thode de regulation selon la revendication 8, appliquee k un 
procede de polymerisation, dans laquelle : 

- on mesure periodiquement I'indice de fluidite en fondu (MI) et/ou la 
masse volumique standard (MVS) du polym&re et/ou sa teneur en 

10 comonomere; 

- Toi^gane de prediction de Talgorithme maitre calcule des predictions 
brutes de MI et MVS en fonction de la temperature au reacteur, des 
concentrations dans le reacteur et des temps de sejour dans les diJferenis 
appareiis du circuit de polymerisation ; 

15 - periodiquement, Vorgane d'adaptation de Talgorithme maitre : 

resynchronise les predictions brutes de MI et MVS, en tenant compte 
du temps ^coule entre la prise des mesures de MI et MVS et 
Tobtention du resultat des mesures, et compare les predictions brutes 
de MI et MVS resynchronisees aux mesures de MI et MVS, 

20 - calcule un paramfetre d'ad^tation multiplicatif kMI applique a la 

prediction brute du MI pour obtenir une prediction calibree du MI, 
et 

calcule un parametre d'adaptation additif kMVS appliqud a la 
prediction brute de MVS pour obtenir une pr&iiction calibree de 
25 MVS. 

15 - Methode de regulation selon Tune des revendications 
prec6dentes, appliqu6e a un procede de polymerisation, dans laquelle on 
6value une ou plusieurs proprietes du polymere en utilisant une technique 
choisie parmi la spectroscopic infrarouge proche (MIR), la spectroscopic 

30 infrarouge a transformee de Fourier (FIIR) et la resonance magnetique 
nucleaire (NMR). 

16 - Methode de regulation selon Tune des revendications 
precedentes, appliquee a un procede de polymerisation, dans laquelle on 
evalue une ou plusieurs proprietes du polymfere en appliquant une relation 
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de correlation preetablie aux resultats de mesures realisees par spectroscopic 
infrarouge proche (NIR) a plusieurs longueurs d'onde predeterminees en 
fonction de la nature du polymere, choisies entre 0,8 et 2,6 mm. 

17 - Procede de synthase d'un produit chimique dans un appareillage 
5 • comprenant au moins un reacteur pouvant etre assimile a un reacteur 

parfaitement melange, regule au moyen de la methode de regulation selon 
Tune des revendications prec&Ientes. 

18 - Dispositif de regulation d*un procede de synthese d'un produit 
chimique dans un appareillage de synthase comprenant au moins un reacteur 

10 pouvant etre assimile a un reacteur parfaitement melange, dans lequel la 
m6thode de regulation est conforme a Tune des revendications 1^16. 
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